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ABSTRACT 


A  method  is  presented  for  the  logical  design  of  single-stage,  combinatorial 
switching  circuits  of  n-variables.  This  method  is  applicable  to  circuits  composed 
of  threshold  devices,  such  as  magnetic  cores,  transistors  with  Kirchoff  adder 
inputs,  parametrons,  etc.  A  study  of  the  constraints  imposed  by  the  form  of  the 
input  portions  of  the  threshold  devices  leads  to  the  definition  of  certain  classes 
of  functions  which  are  physically  realizable  in  a  single  device.  By  the  use  of 
this  method,  arbitrary  switching  functions  of  as  many  as  seven  variables  have 
been  easily  designed  by  hand  computations. 

An  algorithm  for  mechanizing  Boolean  switching  functions,  by  means  of  a 
net  of  magnetic  toroidal  cores,  is  described.  The  algorithm  is  referred  to  as 
Simplex  which,  in  this  application,  is  programmed  for  a  digital  computer. 
Computer-derived  solutions  specify  the  wiring  configuration  for  a  core  or  net 
of  cores  yielding  a  device  for  performing  combinatoried.  logic.  Switching  functions 
are  realizable  in  essentially  one  clock  time.  The  logical  designer  may  synthesize 
a  fimction  directly  from  the  truth  table  without  proceeding  in  the  customary 
manner  of  expressing  the  function  of  Boolean  canonical  form  and  then  attempting 
to  minimize  with  respect  to  hardware  or  other  criteria  by  means  of  algebraic 
manipulation  or  some  mapping  or  charting  technique. 

An  updated  bibliography  is  included. 
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SECTION  I 


ARBITRARY  BOOLEAN  FUNCTIONS  ON  N-VARIABLES  REALIZABLE 
IN  TERMS  OF  THRESHOLD  DEVICES 

INTRODUCTION 

2*^ 

It  is  well  known  that  there  are  2  possible  switching  functions  of  n-variables, 
where  each  variable  can  take  only  the  values  "l"  or  "O. "  Each  truth  function  has 
2^  states.  For  the  purpose  of  designing  arbitrary  functions,  switching  circuits 
for  only  a  small  class  of  these  functions  need  to  be  catalogued,  since  it  can  be 
shown  that  this  class  of  circuits,  under  the  symmetry  operations  of  permutations 

and  complementations  of  the  input  variables,  maps  the  entire  set. 

5 1 

Aiken  lists  all  vacuum  tube  switching  circuits  for  four  variables  in  a  table 
of  402  entries.  Two  other  tables  are  given  which  show  the  transformation  of  the 
function  imder  the  group  of  symmetry  operations  on  the  input  variables,  but  even 
for  four  variables,  the  process  of  transforming  and  classifying  an  arbitrary  function 
is  extremely  tedious,  while  five  variables  require  a  catalog  of  over  one  million 
switching  circuits. 

The  above  classification  is  without  regard  to  the  properties  of  the  physical 
circuit  which  is  used  to  embody  the  function.  Thus,  the  same  "primitive"  set  of 
functions  is  designed  whether  the  circuit  elements  are  vacuum  tubes,  magnetic 
cores,  or  relays,  etc.  If  the  circuit  element  is  a  threshold  device  such  as  a 
magnetic  core,  then  certain  constraints  are  introduced  into  the  problem.  These 
constraints  offer  a  different  class  of  functions,  those  which  can  be  embodied  in  a 

Manuscript  released  by  the  authors  (May  25,  1961)  for  publication  as  an  ASD 
Technical  Report. 
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single  input-composite  of  a  threshold  device.  Such  functions  will  be  called  con¬ 
sistent  with  respect  to  a  threshold  device.  If  only  the  consistent  switching 
circuits  need  to  be  designed  and  allowing  symmetry  operations,  a  table  of  onl’ 

14  entries  is  required  for  all  four  variable  functions,  and  probably  no  more  than 
65  for  all  five  variable  functions. 

The  design  of  an  arbitrary  n-variable  function,  using  threshold  devices  is 
accomplished  in  the  following  manner:  a)  The  truth  function  is  first  partitioned  so 
as  to  separate  the  "one"  bits  into  a  set  of  functions,  each  of  which  is  independently 
consistent,  b)  The  input  composites  for  each  of  these  functions  is  designed,  per¬ 
haps  by  referring  to  a  table  of  consistent  switching  circuits.  c)  Since  the  parti¬ 
tioning  of  the  function  places  each  of  its  "ones"  into  one,  and  only  one,  sub-set, 
the  output  of  the  switching  circuits  need  only  to  be  ORed  together  to  embody  the 
entire  function. 

It  is  the  purpose  of  this  report  to  demonstrate  a  means  of  pa?  dtioning  and 

designing  arbitrary  n-variable  Boolean  functions.  Although  most  of  the  reported 

work  is  with  reference  to  magnetic  cores,  it  is  presented  in  terms  of  the  generalized 

concept  of  a  threshold  switching  device.  For  this  reason  and  because  they  are  ade- 

56 

quately  treated  elsewhere,  the  physical  details  of  magnetic  core  switching  is 
omitted  here. 

This  report  also  covers  a  computer  technique  for  mechanizing  any  arbitrary 
Boolean  function  with  toroidal  cores.  In  other  words,  the  computer-derived 
solution  will  indicate  the  wiring  configuration  for  a  core  or  net  of  cores  yielding 
a  device  for  performing  combinational  logic.  Switching  functions  of  n-variables 
are  realizable  in  one  clock  time. 
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The  technique  is  a  linear  programming  technique  referred  to  as  the  Simplex 
algorithm.  Of  extreme  importance  is  the  fact  that  the  logical  designer  may  syn¬ 
thesize  a  function  directly  from  the  truth  table  without  proceeding  in  the  customary 
manner  of  expressing  the  function  in  Boolean  canonical  form  and  then  attempting 
to  minimize  the  hardware  by  algebraic  manipulation  or  some  other  charting  or 
mapping  technique. 

THRESHOLD  DEVICE  LOGICAL  ELEMENT 

A  threshold  device  can  be  brought  to  either  of  two  distinct  states  depending 
upon  the  magnitude  of  some  physical  quantity  at  its  input.  Although  the  threshold 
of  a  given  device  has  a  distinct  magnitude  (i.e.  in  a  magnetic  core  the  applied  mag¬ 
netizing  force,  Hj^,  must  exceed  the  coercive  force,  H^,  in  order  to  switch)  it  is 
preferable  to  specify  a  range.  That  is  if  the  state  is  not  to  change,  the  strength 
of  the  input  must  be  less  than  a  specified  quantity;  if  it  is  to  change,  the  strength 
of  the  input  must  be  greater  than  another  quantity.  To  use  again  the  example  of  a 
magnetic  core,  if  the  core  is  not  to  switch  from  negative  remanence,  the  applied 
magnetizing  force  H^  ^  0;  or  if  it  is  to  switch  H^  ^  ^c  depends 

upon  the  desired  speed  of  switching.  In  units  of  H^  :  H^  s  0  implies  a  "0,  "  H^  a  1 
implies  a  "1,  "  where  "O"  and  "l"  denote  respectively  the  negative  and  positive 
remanent  states  of  the  core  and  it  is  assumed  that  the  core  always  starts  from 
negative  remanence. 

The  threshold  device  may  be  considered  to  be  composed  of  two  separate 
portions,  an  input  circuit  which  comprises  a  means  of  weighting  and  combining 
the  input  variables,  and  an  output  stage  which  indicates  by  a  definite  change  in 
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state  whether  or  not  the  combined  value  of  the  weighted  variables  exceeds  the 

S6 

threshold.  Following  Karnaugh  °  the  mathematical  description  of  the  input 
circuit  is  called  the  input  composite.  The  input  variables  are  denoted  by 
x^X2  .  .  .  .x^  and  can  take  only  the  values  0  or  1.  All  variables  operate  simul¬ 
taneously  in  the  input  composite.  The  term  x^  is  reserved  for  a  unit  pulse, 
which  together  with  its  weighting  factor,  is  called  the  bias. 

The  effects  of  the  inputs  are  added  algebraically  to  obtain  the  net  effect 
on  the  output  stage  of  the  threshold  device,  and  the  input  composite  may,  there¬ 
fore,  be  represented  in  the  following  form: 


M  =  q  N  X  +  q  N,x,  +  ....  +  q  N  x  (1) 

^o  o  o  1  1  ^n  n  n  '  ' 

Here  q^  =  ±1  and  is  a  positive  integer  or  zero.  In  terms  of  magnetic  cores, 

q^  represents  the  polarity,  the  number  of  turns  of  the  winding,  and  M  the 

net  applied  magnetomotive  force  where  x^  represents  a  constant  current  pulse. 

Now,  since  the  value  of  M  depends  on  the  manner  in  which  the  variables 


x^  are  taken  as  0  and  1,  we  define  a  new  quantity  in  terms  of  a  column  matrix 
whi.  '  will  be  called  the  solution  vector,  S,  the  terms  of  which  are  called  the 
input  values. 


q  N 
o  o 

q  N 
^1-1 


S.  = 

q  N 
n  n 

The  M  above  is  one  term  of  a  vector  which  is  found  by  multiplying  the 
solution  vector  by  the  truth  table  for  n-variables,  treated  as  a  matrix  T .  The 
vector  ^  will  be  called  the  input  composite  vector. 
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X  X  x„  x„ 
0  12  3 


M 


0  0  0 
1  0  0 
1  0 
1  0 


0 
1 

1  0  0 
1 
1 


1  1  1  1  ....  1 


qiNi 

^2^2 


q  N 
n  n 


N 

o  p 

q  N 

o 

+ 

q^Ni 

q  N 
o  o 

+ 

+ 

^2''2 

q  N 

o 

+ 

qiNi 

+ 

‘l2''2 

q  N 
o  o 

+■ 

*13^3 

q  N 
o  o 

+ 

qiN^ 

+ 

‘^2^2 

^3^3 

The  logic  function  ^  is  evaluated  by  applying  the  threshold  condition  to  M . 


SYMMETRY  OPERATIONS 

The  input  circuit  of  a  threshold  device  is  specified  by  its  solution  vector. 
Since  it  is  a  physical  structure,  it  is  invariant  to  complementations  of  and 
permutations  among  the  input  variables.  It  is  easily  shown  that  the  effect  of  the 
symmetry  operations  is  to  permute  the  rows  of  the  truth  table.  Consequently, 
the  terms  of  the  input  composite  vector  are  likewise  permuted  without  changing 
their  values.  One  can,  thus,  postulate  a  2^  x  2*^  permutation-matrix  P 

which  transforms  a  given  M  or  its  corresponding  logic  function,  F,  to  another  F' 
having  the  same  primitive  solution,  or  switching  circuit.  Now,  since  it  is  desir- 
ible  to  avoid  the  use  of  complement  drivers,  and  the  necessity  for  interchanging 
leads,  we  transform  the  solution  matrix,  S,  by  means  of  a  matrix,  which  is 
defined  as  follows: 
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T  • 

S 

=  M 

(2) 

T  • 

c  •  s 

=  =  P-T-S  =  P-M 

(3) 

z  ■ 

c 

=  P  •  T 

(4) 

pre-multiplying 

It  • 

both  sides  by 

T  .  £  =  *  P  •  T 

(5) 

It  can  be  shown  by  partitioning  TT  and  applying  induction  that  ^  is  non- 

singular  and  thus  has  an  inverse,  (T^  •  ^  for  any  number  of  variables,  n. 

C  =  (T^  •  T)"^  •  it  ■  —  ‘ 

Only  certain  P-matrices  are  allowed  because,  although  the  transformation, 

M'  =  ^  is  not  unique,  not  all  of  the  possible  P-matrices  will  lead  to  a 

product,  it  •  i  *  i>  which  is  non-singular.  The  properties  of  the  allowed 
P-matrix  have  not  been  fully  investigated  as  yet.  In  any  case,  we  are  now  in 


a  position  to  describe  ^-matrices  which  will  perform  the  symmetry  operations, 
on  both  the  solution  S  and  the  truth  matrix  by  (4)  above.  An  example  will, 
perhaps,  best  bring  out  the  qualities  of  the  ^-matrix.  Given  a  Boolean  function 
of  3  variables,  F  =  1010110^,  (the  least  significant  bit  is  underlined),  the 
permutation  matrix,  P,  below  transforms  it  to  F'  =  00100111  which  has  the 


same  switching  circuit  and  which  is  listed  in  a  table  of  basic  switching  circuits 


such  as  that  of  Aiken. 

F'  =  P  '  F 


1 

00000001 

1 

00010000 

0 

1 

00000100 

1 

0 

01000000 

• 

1- 

0 

00000010 

0 

1 

00100000 

1 

0 

00001000 

0 

0 

10000000 

1 

(7) 
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The  C  matrix  if  formed  by  the  use  of  formula  (6) 
-C  =  (T^  •  T)'^  •  (T^  •  P  •  ^) 


For  n  =  3,  the  product  •  T  and  its  inverse  are; 

1  -1  -1 
2  0  0 
0  2  0 
0  0  2 

Then  .  P  •  T 


8  4  4  4 
4  4  2  2 
4  2  4  2 
4  2  2  4 


(I*  •  T)  •T 


2 

-1 

-1 

-1 


1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

0 

1 

u 

1 

0 

1 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

“o  0000001 

"l  0  0  0 

00010000 

110  0 

00000100 

10  10 

*8  4  4  4 

01000000 

1110 

2 

4  2  2  0 

00000010 

10  0  1 

4  2  0  2 

00100000 

110  1 

4  02  2 

00001000 

10  11 

1 

10000000 

1111. 

(T  .  T)'^  .  (T,  .  P  •  J)  «  C 
t  ^ 


2  -1  -1  -1 

8  4  4  4 

1111 

-12  0  0 

4  2  2  0 

_ 

0  0  0  -1 

-10  2  0 

4  2  0  2 

0  0-10 

-10  0  2 

4  0  2  2 

0-100 

(8) 


(9) 


Since  post-multiplying  T  by  C  is  the  same  as  pre-multiplying  it  by  P,  let 


us  form  the  product  T  *  ^ 


T'  and  examine  the  resulting  transformation  of 


the  truth  matrix. 
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L- 

0 

0 
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Comparing  ^  and  we  see  that  the  effect  of  the  matrix  ^  has  been  to  inter¬ 
change  columns  and  x^  and  replace  the  variables  by  their  complements  x^’, 
x^',  Xg’.  Indeed,  ^  is  the  product  of  a  permutation  matrix  n  and  a  complemen¬ 
tation  matrix  K. 


TT 

10  0  0 
0  0  0  1. 
0  0  10 
0  10  0 


JK 

1111 
0-100 
0  0-10 
0  0  0  -I 


1111 
0  0  0  -1 
0  0-10 
0-100 


(11) 


Thus,  to  complement  any  variable  x^  without  recourse  to  complement  drivers, 

5  B 

it  is  replaced  by  (x^  -  x^),  this  corresponds  to  Karnaugh's  theorem. 

All  possible  matrices  n  and  K  for  n-variables,  together  with  the  unit  matrix 
form  a  group  which  is  isomorphic  with  the  group  of  isometries  of  an  n-dimensional 
hypercube.  For  n  =  3,  the  group  contains  48  members,  i.  e.,  six  permutations, 
and  eight  complementations. 

We  have,  thus,  shown  that  if  a  permutation  matrix  P  can  be  found  which 
transforms  an  arbitrary  function  to  another  having  a  known  solution  (or  switching 
circuit)  and  if  ]P  has  the  property  such  that  ^  •  IP  •  ^  is  non-singular,  then  a  trans¬ 
formation  C  can  be  found  which  will  convert  the  known  solution  to  the  desired  one. 
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FUNCTION  TRANSFORMATION 


A  function  will  be  called  in  basic  form  when  it  has  been  transformed  by  a 

set  of  symmetry  operations  to  another  function  having  the  same  physical 

switching  circuit,  but  such  that  all  the  true  bits  have  been  moved  as  closely  as 

possible  to  the  least  significant  bit  end  of  the  function.  The  required  set  of 

symmetry  operations  can  be  selected  in  the  following  manner. 

Consider  a  truth  function  on  n-variables.  The  nth  variable  is  false  for  the 

first  half  and  true  for  the  second  half.  Count  the  number  of  true  states  of  the 

function  in  the  second  half  of  the  function  and  subtract  from  the  number  of  true 

bits  in  the  first  half  of  the  function.  Call  this  number  A  .  If  A  <0,  there  are 

n  n 

more  true  states  in  the  second  half  of  the  function  than  in  the  first  half.  Thus, 

complementing  the  variable  is  one  of  the  required  symmetry  operations.  A 

set  of  numbers  A^,  1  i  s:  n,  is  found  by  subtracting  the  number  of  states  of 

the  function  which  are  true  when  the  variable,  X.  is  true  from  the  number  of 

1 

states  which  are  true  when  the  variable  is  false.  With  this  set  of  numbers,  it 
is  possible  to  describe  the  required  set  of  symmetry  operation. 

Consider  now  a  function  in  basic  form  which  has  k  true  states  followed  by 
(m  -  k)  false  states  where  m  =  2*^.  The  following  statements  concerning  the  A 
and  solution  vectors  can  be  inferred  from  the  form  of  the  function.  Proofs  are 
omitted  here. 

1)  Parity  is  conserved  for  all  k  among  k  and  A..  This  is  also 
true  for  any  function. 
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2)  The  can  have  the  following  values  only: 

0  i  Aj  ^  1 

0  ^  A2  ^  2 

0  ^  Ag  i  4 

0  ^  A,  8 

4 

0  ^  A  <  2””  ^ 
n 

3)  ij.i  - 

4)  All  k  true  out  of  m  functions  with  no  gaps  are  realizable. 

5)  The  minimal  solution  vector  will  have  the  following  form: 

qo  =  +1 

q.  =  -  1  for  1  i  s  n 

1 

N,  s  N„  .  .  .  .  ^  N  s  N 
12  no 

6)  The  statements  are  equivalent  that  the  given  function  is  in  basic 
form  when  a)  the  Ai  are  in  monotonically  increasing  order  or  b) 
the  solution  if  it  exists  is  as  defined  in  5  above. 

We  shall  now  state  by  analogy  with  the  function  studied  that  any  arbitrary 
function  is  in  basic  form  when  its  A^  are  all  positive  and  in  monotonically  in¬ 
creasing  order  as  i  increases.  The  set  of  symmetry  operations  required  to 
transform  the  arbitrary  function  to  basic  form  is  identical  with  the  operations 
required  to  transform  the  set  A^^  to  basic  form.  Since  complementations  and 
permutations  do  not  in  general  commute,  all  variables  with  negative  A^  are 
first  complemented  and  then  the  indicated  permutations  are  performed. 
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With  some  practice,  the  transformation  of  the  functions  for  any  symmetry 
operation  becomes  obvious.  For  example,  interchanging  and  permutes 
the  terms  of  the  pairs;  1,2;  5,6;  9,10;  13,14;  etc.  of  the  truth  table.  Com¬ 
plementing  x^  permutes  the  terms  of  the  pairs:  0,  1;  2,3;  4,5;  etc. 

The  C  matrix  required  to  transform  the  solution  in  basic  form,  to  the 
required  solution,  S,  is  found  as  follows.  If  ^  is  the  vector  for  the  given  function 

T 

and  ^  is  the  vector  for  the  same  function  in  basic  form  one  can  write  a  matrix 
D  such  that 

D  •  a'  =  ^ 

Since  the  solution  has  the  order  r.+ 1,  the  D  matrix  is  converted  to  the  C 
matrix  by  bordering  with  a  row  and  a  column.  For  example: 

Example  1 

F  =  00101010;  A^=-3  A2=l  A^^l 

F'  -  0  0  0  0  0  1  1  ^  ;  a'j  =  1  a'^  =  1  A'^  *  3 

D  .  a' 

T  0  -T[  [i 

0  1  0|  1 

10  0,  3 

r  "* 

I  1  0  0  1 

'•0  0  0  -1 
—  "  0  0  10 
0  1  0  oj 

To  perform  complementation,  every  column  which  contains  a  negative 
entry  must  have  a  one  in  the  added  row. 


A 
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» 

I 


s 


c_  .  S' 

10  0  1  2 

000-1  -1 

0  0  10  •  -1 
0  10  0  -2 


It  is  apparently  true  that  the  set,  Aj,  uniquely  specifies  the  function  if  it  is 
realizable  in  a  single  input  composite,  and  it  is  obvious  from  the  way  they  are 
defined  that  the  absolute  values  of  the  members  of  the  set  for  a  given  function 
are  invariant  to  symmetry  operations  on  the  function,  although  their  order  may 


be  interchanged.  If  the  function  is  not  realizable  in  a  single  input  composite. 


the  set,  A^,  does  not  uniquely  specify  it,  for  example  see  figure  la. 


A  B 
1  1 
1  1 
1  1 
1  1 


10  0 
10  0 
10  0 
1  0  0 


1  0  1  0  0  0  0 
0  1  0  0  0  1  0 
0  1  0  0  0  0  1 
1  0  0  1  0  0  0 


11  10  10  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 


0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 


Figure  la. 

Non  Realizable  Functions 


Figure  lb. 

Partitioned  Functions 
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The  non-realizable  fimctions  A  and  B  in  figure  la  have  the  same  set,  A^, 
but  neither  can  be  transformed  into  the  other  by  any  set  of  symmetry  operations. 


showing  that  the  A^  does  not  specify  uniquely  a  non-realizable  function.  Partitioning 

the  functions  A  and  B  in  figure  lb  it  is  seen  that  summing  over  the  like  terms  of 

^  vectors  of  the  realizable  functions  making  up  the  partitions  produces  the  A 

vector  of  the  non-realizable  function. 

Given  a  realizable  function,  if  A^  =  A^  permuting  the  variables  and 

will  cause  certain  true  states  of  the  function  to  commute  without  changing  the 

function,  since  it  implies  that  q  N  *  q,  N,  .  If  the  function  is  not  realizable, 

a  a  b  b 

however,  permuting  the  inputs  X^  and  Xj^  may  change  the  function  even  if  A^  =  A^^. 
Thus,  if  the  function  is  not  realizable,  there  is  no  guarantee  that  ordering  its 
A^  monotonically  will  cause  the  true  states  to  be  clustered  most  densely.  Never¬ 
theless,  the  function  will  be  considered  in  basic  form  when  its  set  of  A^  is 
monotonically  ordered  w.  r.  t.  increasing  i. 


CONSISTENT  FUNCTIONS 

If  we  examine  a  truth  table  for  n-variables,  we  see  that  each  span  of  four 
rows  repeats  the  truth  values  for  the  variables  x  and  x„.  The  variable  x„ 

X  «  o 

becomes  true  in  the  second  span  and  remains  true  for  the  entire  span,  likewise 
with  x^,  in  the  third  span.  In  fact,  each  variable,  x^,  3  s  i  s  n,  must  be 
constant  in  any  arbitrary  span  of  fourxows.  Thus,  any  variations  in  the  values 
of  the  terms  of  the  input  composite  vector  within  any  span  of  four  terms  must  be 
due  to  the  variables  x^  and  x^.  Within  any  span,  the  presence  of  true  variables 
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whose  indices  are  3  or  greater  only  affect  the  net  bias  and  establish  an  average 
level  for  that  span  as  can  be  seen  in  figure  2.  Consequently,  it  is  ii  Jtructive 
to  study  the  properties  of  the  set  of  functions  of  two  variables. 

Of  the  16  functions  of  two  variables,  14  of  them  are  consistent  as  defined  in 
the  introduction.  With  a  little  reflection,  it  is  possible  to  write  unique,  minimal 
solutions  for  each  of  them.  These  functions  are  shown  in  Table  I,  where  the 
decimal  heading  indicates  the  function,  i.e.  one  is  0001,  five  is  010^^,  etc.  Six 
and  nine,  the  two  alternating  symmetric  functions  (EXCLUSIVE  OR  and  its 
inverse),  cannot  be  embodied  in  a  single  input  composite  and  are,  therefore, 
inconsistent.  It  is  easy  to  see  that  this  is  true  by  referring  to  figure  2.  That 
is,  for  function  "six,  "  if  both  q.N,  +  q  N  ^  1,  and  q  N  +  q„N_  ^  1  where 
q  N  ^0,  then  surely  q  N  +  q.N  +  q„N„  ^  1,  and  similarly  for  "nine." 

OO  O  O  1.  X  Z  (6 

As  an  example  of  the  manner  in  which  solutions  of  'dgher  order  functions 

can  be  built  up  from  the  two-variable  solutions,  compare  the  solution  for  function 

"seven"  with  the  solution  for  function  "one.  "  We  see  that  the  two  are  identical 

if  the  x^  input  value  of  the  former  is  reduced  by  one.  If  we  take  the  three -variable 

function  f  =  0001,  0111,  the  first  span  can  be  represented  by  function  "seven" 
o 

and  the  second  by  function  "one.  "  Referring  to  figure  2,  in  the  second  spEui, 
the  net  bias  term,  which  is  constant  throughout  the  span  is  q  N  +  q^N-.  If 
we,  therefore,  make  the  input  value  of  x  ,  q,N„  “  -  1,  then  the  solution  for  the 
first  span  fits  the  second  span  and  the  input  composite  is  M  =  2x^  -  Ix^  -  Ix^  -  Ix^. 
Plotting  the  successive  values  of  M  (the  terms  of  the  input  composite  vector)  we 
get  the  input  composite  profile  which  aids  in  visualizing  the  described  process. 

(See  figure  3)  Thus,  the  profiles  of  the  minimal  solutions  of  the  functions  "seven" 
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■ 

1 

1 

0 

1 

1  . . 

..  1 

(qN+qN+q  ,N 
o  o  n  n  n-1  n-1 

1 

0 

1 

1 

1  . . 

. .  1 

(qN+qN+q  .N  ,.. 
o  o  n  n  n- 1  n- 1 

13«3>  *<'2'‘2 

1 

1 

1 

1 

1  . . 

. .  1 

(qN+qN+q  ,N 
o  o  n  n  n-1  n-1 

‘>3''3’*‘>l”l*‘‘2''2 

FIGURE  2  -  A  Portion  of  a  Truth  Table  for 

N-Variables 
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No. 

0 

1 

2 

3 

4 

5 

7 

8 

10 

11 

12 

13 

14 

15 

0 

1 

0 

1 

0 

1 

1 

0 

0 

1 

0 

1 

0 

1 

Fn. 

0 

0 

1 

1 

0 

0 

1 

0 

1 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

0 

1 

0 

1 

2 

-1 

0 

1 

0 

1 

0 

1 

S. 

0 

-1 

1 

0 

-1 

-1 

-1 

1 

1 

1 

0 

-1 

1 

0 

0 

-1 

-1 

-1 

1 

0 

-1 

1 

0 

-1 

1 

1 

1 

0 

TABLE  I.  Consistent  Solutions  of  Two  Variable  Functions 

and  "one"  are  identical  in  form,  but  that  the  latter  is  shifted  one  unit  to  the  left. 
Likewise,  the  two  spans  of  the  function  F  are  identical  in  form  with  the  lower 
span  shifted  one  unit  to  the  left  with  respect  to  the  upper  span,  in  accordance 
with  the  input  value  of  x^.  The  function  F  is,  therefore,  consistent. 

Functions  such  as  "seven"  and  "one"  will  be  called  compatible  with  one 
another  if  both  functions  can  be  derived  from  the  same  solution  merely  by  a 
change  in  the  bias.  It  can  be  shown  that  other  compatibilities  exist  if  certain 
equivalent  modified  two-variable  solutions  are  admitted.  The  only  allowed 
modiiications  of  the  minimal  solutions  are  in  the  terms  N  and  N_  of  the  solution 

1  M 

vector.  The  signs  of  the  termsj,  are  not  allowed  to  change.  The  can  only 
be  increased,  since  we  start  with  a  minimal  solution  and  reducing  the  input 
value  of  either  variable  to  zero  reduces  the  order  of  the  solution.  Likewise, 
if  the  input-value  of  x^  is  zero,  we  cannot  alter  it  without  more  information. 
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FIGURE  3  -  Function  Profiles 


2 

-1 

-1 

-1 


since  the  value  of  is  not  known.  We  now  investigate  whether  a  modified,  or 
non-minimal  solution,  which  may  represent  a  different  function  can  be  restored 
to  representing  the  original  function  merely  by  a  change  in  bias. 

If  a  modified  two-variable  solution  is  formed  from  a  minimal  two-variable 
solution  by  increasing  the  input  value  of  one  of  the  two  variables,  it  can  be 
seen  by  the  first  span  of  figure  2  that  the  value  of  the  input  composite  for  either 
pair  of  terms  M.  and  M„  or  M„  and  M_  will  be  altered,  without  altering  the 

i.  O  A 

difference  between  them  because  the  same  amount  is  added  to  each.  Thus,  if 
the  truth  values  of  either  or  both  members  of  one  pair  of  terms  change,  their 
truth  values  can  be  restored  by  a  compensating  change  in  the  bias.  This  compen¬ 
sating  change  will  not  alter  the  truth  values  of  the  unaffected  terms  for  the 
following  reasons:  Consider  first  the  terms  and  Mg. 

(a)  If  both  terms  are  greater  than  one,  and  if  the  input  value  of  Xj^  or  Xg 


-17- 


is  altered  so  as  to  cause  either  Mj^  or  Mg  to  become  zero  or  less  (the  of  the 
variable  in  question  must  be  -  1),  then  in  "restoring"  by  a  bias  change,  which 
must  be  in  the  positive  direction,  the  previously  unaffected  term  is  made  more 
positive  so  its  truth  value  is  unchanged. 

(b)  If  both  terms  were  zero  or  less,  by  the  same  reasoning,  the  truth  value 
of  the  unaltered  term  is  not  affected  by  restoring. 

(c)  If  one  term  is  one  or  greater  and  the  other  is  zero  or  less,  the  truth 
values  of  terms  and  Mg  cannot  be  altered  by  altei  ing  or  Ng  since  the  signs 
q  and  q  are  opposite.  Note  that  if  either  N  or  N„  is  zero,  it  cannot  be  altered 

i.  A  1  A 

since  we  do  not  know  the  sign  of  its  q^. 

Now  consider  the  term  M^,  the  truth  value  of  which  is  most  likely  to  change 
with  restoring  by  bias  compensation  since  it  depends  only  on  N^: 

(a)  If  M  s  0  and  the  solution  is  to  be  altered  so  as  to  cause  either  M,  or 

o  1 

Mg  to  exceed  the  threshold,  it  can  only  be  done  if  the  appropriate  input  value  was 

initially  positive.  Restoring  will,  therefore,  only  make  M^  more  negative  and 

will  not  alter  its  truth  value,  f  . 

o 

(b)  If  M  SO  cind  either  M,  or  M„  is  to  be, made  less  than  zero,  it  cannot  be 

O  i.  ^ 

done  by  the  allowed  alterations,  since  the  appropriate  q^  must  be  initially  positive. 

(c)  If  M^  2  1  the  same  reasoning  shows  that  the  only  possible  alteration  in 

the  term  Mj^  or  Mg  is  such  as  to  cause  its  value  to  go  from  one  or  greater  to  zero 

or  less,  changing  its  truth  value  from  a  "l"  to  a  "O.  "  Therefore,  restoring  the 

truth  value  of  the  term  M,  or  M„  will  not  alter  the  truth  value  of  M  . 

12  o 

Thus,  it  has  been  shown  that  a  minimal  two-variable  solution  which  represents 
one  two-variable  function  can  be  altered  so  as  to  represent  a  different,  though 
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related,  function  by  altering  the  input  value  of  one  or  both  of  the  variables, 

X.,  x„,  or  the  bias  x  and  then  "restored"  to  representing  the  original  function 
i.  o 

solely  by  a  change  in  the  bias  input  value  These  related  functions,  sharing 

a  non-minimal  solution  are  also  compatible. 

As  an  example,  let  us  solve  the  three-variable  function^  =  0001010_1. 
Dividing  it  into  two  four-bit  spans,  noting  their  minimal,  two-variable  solutions 
from  Table  I,  and  generating  their  profiles,  we  have  figure  4. 
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FIGURE  4  -  Sum  of  Minimal  Profiles 


It  is  readily  seen  that  no  possible  alteration  of  the  input  values  of  the  biases 
of  either  two-variable  solutions  in  figure  4  can  make  either  function  fit  the  other 
function.  Now,  if  we  add  like  terms  of  the  solution  vectors,  or  add  across  the 
two  profiles,  we  obtain  the  non-minimal  solution  and  profile  shown  in  the  third 
column  of  figure  4,  which  still  represents  the  upper  span  function.  The  fourth 
column  shows  that  by  adding  -  1  to  the  bias,  the  same  solution  and  profile  can 
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also  represent  the  second  span  of  the  function  F.  This  indicates,  as  in  the 


previous  example,  that  the  input  value  of  variable  x^,  =  -  1  and  the  solution 

for  the  function  F  is  as  shown  in  figure  5.  If  the  two  corresponding  terms  of  two 
minimal  solutions  have  opposite  signs,  the  order  of  the  solution  is  reduced  and 
the  sum  cannot  represent  either.  Such  solutions  are  not  compatible. 


T  S  =  ^  Profile  F 

lOOo]  2  [  2!  _ 1  1 

1100  -2  0  0 

1010  -1  1  !  I  1 

1110  -ij  -1  I  i  0 

1001  1  !  I  1 

1101  -1  Llj  0 

1011  0  I  0 

_llllj  -2  1  ^  0 


-2-10+1+2 

FIGURE  5  -  A  consistent  Three-Variable  Function 

and  Solution 


A  A'  B  B'  C  CD  D' 


FIGURE  6  -  Consistent  Minimal  Profiles 


The  fourteen  consistent,  two-variable  functions  may  now  be  tabulated  and 
classified  with  respect  to  their  compatibility  with  one  another.  A  convenient 
way  to  do  so  is  on  the  basis  of  their  profiles,  figure  6. 
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Figure  7  -  CIRCLE  OF  COMPATIBILITY 

EACH  SOLUTION  IS  COMPATIBLE  ONLY 
WITH  ITS  NEAREST  NEIGHBORS. 

The  function  is  written  between  the  values  M  =  0  and  M  =  1.  Where  the 
same  profile  can  represent  two  functions  by  shifting,  they  are  displayed  accord¬ 
ingly.  The  letters  are  arbitrarily  assigned  for  convenient  reference.  The 
compatibilities  of  the  minimal  functions  are  shown  in  figure  7.  The  empty  set 
is  called  "0,  "  the  filled  set  I  and  they  are  compatible  with  all  other  functions. 
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NON-CONSISTENT  FUNCTIONS 


It  will  be  recalled  that  the  functions  "six"  and  "nine"  the  EXCLUSIVE-OR 
and  its  inverse  are  not  consistent.  A  two-variable  function  can  be  thought  of 
as  being  composed  of  four  functions  of  (n-2)  =  0  variables,  and  the  profile  of  a 
zero-variable  function  is  obviously  a  straight  line  at  zero,  or  at  one.  Let  us 
construct  a  hypothetical  profile  for  the  function  "six"  by  shifting  the  zero  variable 
profile  to  accommodate  its  four  quarters,  see  figure  8. 


First  quarter 

1  0 

-  shift  =  +1 

Second  quarter 

1  1 

Third  quarter 

.  1  1 

shift  =  -  1 

Fourth  quarter 

0 

Figure  8  -  Hypothetical  EXCLUSIVE-OR  Profile 

We  see  that  to  accommodate  the  second  quarter  it  is  required  to  shift  the 
(n-2)-variable  profile  by  S  =  +  1  unit,  and  to  accommodate  the  fourth  quarter, 
the  profile  must  be  shifted  by  S  =  -  1  unit.  Thus,  the  absolute  value  of  the 
difference  in  shift  between  the  quarters  in  the  first  half  and  the  quarters  in  the 
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second  half,  6^  =  2.  Referring  to  figure  6,  we  see  that  for  all  two-variable 
consistent  minimal  profiles,  the  difference  in  shift  6j,  =  0.  Referring  now  to 
the  profiles  in  the  first  two  columns  of  figure  4,  the  first  profile  represents  the 
upper  half  and  the  second  the  lower  half  of  the  profile  of  the  3-variable  function, 
0001010  The  shifts  between  the  quarters  for  the  profiles  shown  are 
zero  for  the  first  half  and  -  1  for  the  second  half,  so  that  the  difference  in  shifts 
6g  =  I  0  -1  I  =  1,  this  function  is  consistent,  and  can  be  realized  in  a  single 
input  composite  as  shown  in  figure  5.  It  appears  to  be  a  general  criterion  that 
if  the  absolute  value  of  the  difference  of  shift,  6  S:  2,  the  function  is  not  con- 
sistent  and  requires  more  than  one  input  composite  for  its  realization.  This 
criterion  will  be  discussed  further  under  partitioning. 

THE  PARTITIONING  OF  FUNCTIONS 

Having  examined  some  of  the  properties  of  the  minimal,  consistent  two- 
variable  solutions,  we  are  now  in  a  position  to  consider  a  method  of  partitioning 
the  true  terms  of  an  arbitrary  n-variable  truth  function  into  a  set  of  functions 
each  of  which  is  consistent,  in  accordance  w’th  the  design  procedure  laid  out  in 
the  introduction.  It  was  shown  in  the  discussion  of  consistent  functions  that  any 
term  to  term  variations  of  the  M  vector  within  a  span  of  four  terms  must  be  due 
only  to  variables  and  x^.  Thus,  if  the  function  is  divided  into  spans  of  four 
terms  as  in  figure  2,  the  minimal  two-variable  solutions  for  each  span  considered 
independently  must  be  compatible  with  alt  others  if  the  function  is  consistent. 

Note  that  compatibility  is  not  transitive,  i.  e.  the  fact  that  B  is  compatible  with 
A  and  A  with  C  does  not  imply  that  B  is  compatible  with  C. 

-23- 


It  was  conjectured,  at  first,  that  permuting  every  possible  pair  of  variables 

with  the  variable  x  and  x  successively,  leading  to  n(n-l)/2  related  functions, 

1  « 

(that  is,  related  by  having  the  same  assumed  physical  embodiment,  but  altered 
by  symmetry  operations),  and  checking  each  of  these  derived  functions  for 
compatibility  among  its  four  bit  spans  would  be  a  sufficient  test  for  realizability. 
This  conjecture  turns  out  to  be  false  for  more  than  five  variables.  Since  the 
condition  is  mathematically  necessary,  however,  it  can  be  used  to  delete  incom¬ 
patible  bits  from  a  non-realizable  function.  Although  if  the  resulting  partition 
function  has  more  than  five  variables,  it  may  still  not  be  realizable.  This 
process  is  called  partitioning  by  the  method  of  sieving,  and  a  simple  example 
is  given  in  Table  II. 

In  the  light  of  the  proposed  criterion  for  realizability,  it  is  possible  to 
determine  why  the  2-variable  sieving  test  breaks  down.  The  criterion  essentially 


states  that  if  the  input  value  of  at  least  one  of  the  variables,  x^,  is  not  independent 

of  the  rows  of  the  truth  table  but  requires  in  one  portion  and  in  another 

portion  q^^  N^2'  ^  absolute  value  of  the  difference  between  the  input 

values,  6_  =  I  q.,  N.,  -  q.„N.„  I  i  2,  the  function  is  not  realizable.  The 
S  '  ^]1  jl  ^j2  j2  ' 


difference,  6g,  is  determined  by  the  shift  required  to  accommodate  a  suitably 
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II 


I  FIVE  VARIABLE  FUNCTION 


First  Partition 


10  1110  0  1 
l-Fl  1  0  0  0  0 
1  -F  1  0  0  O-F  0 
1  0  0  0 -Ir -t -Ir  0 
I  BCACC'A'C 
B'  C 


First  Sieve 


Solutions 


The  checked  columns  have  compatible 
solutions.  The  true  bits  of  the  EXCLU- 
SIVE-OR  functions  in  columns  2  and  5 
are  treated  as  separate  functions.  The 
lower  bit  of  column  5,  representing  a  C 
solution  is  stricken  out. 


Perm. 


0000  1001 
10  10  10  11 
0000  0000  Second  Sieve 
1  0  0  0  1  0  1  1 
D'OB  0  B  OD'B  Solutions 


All  columns  compatible,  all  bits  pass 


second  sieve. 

0000  4-04'-!' 
0  0  10  10  10 
0000  0000 
10  10  10  11 
C  OD'O  B  OB  C 
C 

C'OD'O  D'OL’C 


Solutions 


Perm. 


Striking  out  the  bits  in  the  first  row 
converts  the  solutions  to  those  shown 
below  the  line.  All  of  which  are  com¬ 
patible  . 


TABLE  II 
(Cont'd) 

PARTITIONING  OF  A  FIVE  VARIABLE  FUNCTION 

Taking  the  columns  of  the  third  sieve  in  order  from  the  right,  the  first  partition 

results  in  the  function  marked  A,  which  is  consistent.  Its  solution  is: 

-3 

3 

^  =  2 
-1 
-1 

The  successive  partitionings  proceed  in  a  similar  manner. 


chosen  common  profile  to  different  portions  of  the  function. 

There  are  at  least  three  ways  in  which  6g  ^  2  can  occur: 

(1)  qji  =  -1.  qj2  =  +1 

(2)  q.^  N.^  <  W<  qja  =  ‘Ik  =  ^ 

‘Ik  ^k  ^  ‘Ijl  Njl  ^  ‘»J2  <Nj2  +  2)  ^  q^  or 

‘Ik  ^k  ^  ‘Ijl  Nji  ^  qj2  <Nj2-2)  ^  q^N^ 

for  q,  =  q.,  =  q.„  =  q,  =  ±1  and  such  that  no  other  input  has  a  value  intermediate 

K  Ji.  -C 

between  those  of  x,  and  x , . 

k  I 

Condition  (1)  is  equivalent  to  complementing  one  of  the  inputs  of  a  minimal 
compatible,  2-variable  function.  For  example,  complementing  the  variable 
in  a  C  profile  produces  the  B  profile.  This  change  is  detected  in  at  least  one  of 
the  sieves. 
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Condition  (2)  in  the  sieves  is  equivalent  to  a  permutation  of  the  variables 

X.  and  X  .  Thus,  permuting  the  variables  X,  and  X  in  A-C  profile  gives  rise 
3  1  2 

to  a  C-D  profile  and  since  compatibility  is  not  transitive,  condition  (2)  is  also 
detected  in  the  sieving  process. 

Condition  (3)  cannot  be  detected  in  a  sieve  since  only  relationships  such  as 

q.  N  >  q.  N.,  q.  N.  >  q  N  ,  q  N  >  q  N  would  result,  and  would  indicate 
kk  13  13  •C'Ckk 

an  apparent  compatibility. 

The  two  variable  sieving  method  is  equivalent  to  the  two  monotonicity  test 
77  95 

for  realizability  '  .  Similar,  tests  could  be  described  which  would  test  higher 

ordered  monotonicities  for  functions  of  more  than  five  variables  but,  since  it 
is  known  that  complete  monotonicity  is  not  sufficient  for  realizability  and  because 
the  labor  involved  becomes  great,  it  is  not  considered  worthwhile. 

A  partitioning  method  can  be  based  on  the  profile  shift  criterion  mentioned 

above  by  taking  successively  higher  ordered  profiles  as  shown  in  figures  9a  and 

* 

9b.  This  method  partitioned  Moore's  function  into  two  realizable  functions, 
the  one  having  only  a  single  true  state.  The  process  also  resulted  in  solutions 
for  each  of  the  functions. 

CLASSIFICATION  OF  REALIZABLE  FUNCTIONS 

All  of  the  four-variable  single-input-composite  switching  circuits  or  solutions 
are  tabulated  in  Table  III  along  with  the  function  in  basic  form  which  it  represents. 

Moore's  function  is  a  completely  monotonic  function  of  12  variables  which 
is  not  realizable  in  a  single  threshold  device. 
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5  Variable  .  .  Minimal  Common  Minimal  Common  Minimal  Common  3  Variable  Profile 

Function  2  Variable  Profile  3  Variable  Profile  After  Elimination  of  Bit  #18 


New  5  Variable  Function  Minimal  Common  Solution  Function  Solution 

After  Elimination  of  Bit  #18  '  ^  4  Variable  Profile  of  A  B  of  B 


29 


Figure  9b  -  Profile  Shift  Criterion  for  Realizability 


k  n 

The  symbol  )  stands  for  k  true  out  of  m  =  2  states  of  the  truth  table.  The 

solutions  were  extracted  from  a  table  of  all  four  variable,  combinational,  linear- 

input  switching  circuits  derived  by  R.  C.  Minnick  by  means  of  the  simplex 

algorithm  of  linear  programming  using  the  Burroughs  220  computer. 

Admitting  the  operation  of  inversion  along  with  the  symmetry  operations,  the 

fourteen  functions  tabulated  suffice  to  map  every  realizable  function.  Inversion 

is  equivalent  to  subtracting  the  given  function  from  one  which  is  true  for  every 

state.  Whenever  the  number  of  true  states  of  the  realizable  function  is  greater 
n*  1 

than  2  ,  it  is  advantageous  to  invert.  The  solution  is  inverted  by  a  vector 

subtraction,  for  example: 


1 

- 

F 

s 

F 

1 

-  1 

^2 

0 

2 

-2 

0 

-3 

3 

0 

- 

1 

S' 

-1 

0 

-2 

2 

0 

1 

-1 

- - 

Using  the  table  for  four  variable  functions,  it  is  possible  to  build  up  a  table 

0 

for  five  variables  by  attaching  one  function  to  the  end  of  another.  A  (32)  five 

1  5 

variable  basic  function  is  formed  by  attaching  a  (15)  function  to  the  end  of  a  (5^0) 

function.  There  are  two  choices  of  (j^)  functions  and  we  must  take  the  one 

5 

having  a  true  state  in  its  lower  half.  That  is,  the  one  on  the  right  of  the  (^0) 
column  of  Table  IV.  Remembering  that  the  input  value  of  the  new  variable, 
is  negative  and  that  i  if  the  5 -variable  function  is  to  be  in  basic  form, 
the  first  half  of  the  lower  4-variable  function  must  be  contained  in  the  second 
half  of  the  upper  function.  Thus,  since  there  are  no  true  states  in  the  lower  half 
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5 

of  the  {]^0)  function  on  the  left  hand,  no  true  states  can  be  added.  Furthermore, 

since  the  same  4-variable  profile  must  generate  both  halves  of  the  5-variable 

function,  only  compatible  changes  are  allowed.  These  criteria  are  sufficient  for 

all  five  variable  basic  realizable  functions,  but  not  for  six  variables,  and  in 

general,  it  is  not  known  how  to  generate  every  realizable  function  for  more  than 

six  variables.  Table  IV  lists  all  realizable  functions  of  5 -variables  in  basic  form 

along  with  their  A  and  solution  vectors.  Note  that  all  realizable  functions  of 

fewer  than  5  variables  are  automatically  included. 

The  A  vector  indicates  the  number  of  essential  variables  in  the  function; 

if  any  of  the  terms  =  k  where  k  is  the  number  of  true  states  of  the  function, 

it  means  that  the  corresponding  variables  in  the  basic  solution  are  present 

only  to  suppress  true  states.  Thus,  they  are  given  input  weights,  *  ”  ^o‘ 

A  function  of  five  or  fewer  essential  variables  can  be  tested  for  realizability 

by  reference  to  the  table.  First,  determine  the  A  vector  for  the  given  function. 

If  the  function  has  five  or  fewer  essential  variables,  the  table  can  be  used. 

Look  for  a  matching  set  of  A^  in  the  table,  ignoring  negative  signs  and  the 

Ic  1 

order  of  the  terms,  under  the  columns  headed  by  (32).  If  k  >  2  look  in  the 

k*  itc  n 

columns  headed  (32^'  where  k  =  (2  -  k)  for  the  basic  form  of  the  inverted 

function.  The  solution  vector  of  the  basic  function  (if  the  given  function  is 
realizable)  is  converted  to  the  desired  solution  as  in  example  1.  But,  if  the 
function  was  inverted,  the  basic  solution  is  first  inverted  and  then  each  input  is 
complemented  before  deriving  the  ^  and  matrices. 
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Primitive  Functions  of  up  to  Rve  Variables 
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Finally,  the  number  of  entries  in  Tables  III  and  IV  can  be  reduced  by  applying 
the  concept  of  profile  shifting.  Suppose  that  the  bias  input  value  is  not  constant, 
but  a  quantity  that  can  be  varied  positively  and  negatively  in  integer  steps.  The 
function  that  is  generated  will  then  depend  upon  the  bias  setting,  and  the  number 
of  different  functions  which  can  be  generated  depends  upon  the  extent  and  com¬ 
plexity  of  the  profile.  It  can  be  shown  that  by  allowing  the  operation  of  profile 
shifting,  only  two  special,  non-minimal  designs  are  required  for  all  realizable, 
basic-form  functions  of  three  variables,  three  designs  for  four  variables,  and 
eight  for  five  variables. 


-34- 


SECTION  n 


THE  USE  OF  THE  SIMPLEX  ALGORITHM  IN  THE  MECHANIZATION 
OF  BOOLEAN  SWITCHING  FUNCTIONS  BY  MEANS  OF  MAGNETIC  CORES 

INTRODUCTION 

In  recent  years  many  magnetic  devices  have  been  introduced  which  are  used 

to  mechanize  switching  functions.  One  class  of  devices  which  has  received  wide 

attention  is  the  multi-aperture,  or  multi-path,  structures,  such  as  the  Trans- 
19  42 

fluxor,  the  MAD's,  and  the  Laddie,  Their  common  property,  in  addition  to 
using  magnetic  material  with  a  rectangular  hysteresis  loop,  is  that  some  or  all 
of  the  logical  functions  are  performed  by  controlling  the  actual  switching  path 
through  the  structure.  The  philosophy  is  to  achieve  as  much  logic  as  possible 
within  the  structure,  thereby  reducing  component  count  and  interwiring. 

The  question  arises,  however,  as  to  whether  the  simple  toroidal  core  has 

been  given  sufficient  consideration.  A  study  of  Karnaugh's  paper,  "Pulse 

56  65 

Switching  Circuits  Using  Magnetic  Core,  "  and  an  internal  report  by  R.  C.  Minnick 

revealed  extremely  interesting  logical  properties  of  the  simple  toroidal  core 

based  on  the  exploitation  of  its  magnetic  threshold.  For  example,  the  logical 

threshold  function  can  be  mechanized  in  a  single  core,  including,  of  course,  the 

AND  and  the  OR  functions  which  are  special  cases  of  the  threshold  function.  This 

leads  directly  to  the  consideration  as  to  how,  if  possible,  does  one  synthesize 

any  arbitrary  Boolean  function  with  a  single  core  or  perhaps  several  cores. 
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This  section  reports  on  a  computer  technique  for  mechanizing  any  arbitrary 


Boolean  function  with  toroidal  cores.  In  other  words,  the  computer-derived 
solution  will  indicate  the  wiring  configuration  for  a  core  or  net  of  cores  yielding 
a  device  for  performing  combinational  logic.  Switching  functions  of  n  variables 
are  realizable  in  essentially  one  clock  time. 

The  technique  is  a  linear  programming  technique  referred  to  as  the  Simplex 
algorithm.  Of  extreme  importance  is  the  fact  that  the  logical  designer  may 
synthesize  a  function  directly  from  the  truth  table  without  proceeding  in  the 
customary  manner  of  expressing  the  function  in  Boolean  cannonical  form  and 
then  attempting  to  minimize  the  hardware  by  algebraic  manipulation  or  some 
other  charting  or  mapping  technique. 

It  will  be  interesting  to  note  that  this  technique  is  not  limited  to  magnetic 
cores  but  should  apply  to  any  device  with  a  physical  threshold. 

NOTATION 

The  schematic  representation  of  magnetic  circuits  is  facilitated  by  the  use 
of  mirror  symbols.  In  this  paper  cores  (toroids)  are  represented  by  heavy 
vertical  line  segments,  wire  leads  by  horizontal  line  segments,  and  windings  by 
45°  mirror  symbols  at  the  intersections  of  the  cores  and  leads,  figure  10.  The 
sense  of  the  magnetic  field  associated  with  a  current  in  a  given  winding  is  ob¬ 
tained  by  "reflecting"  the  current  in  the  winding  mirror  symbol.  For  example, 
if  the  core  is  initially  in  the  reset,  application  of  current  1^^  to  winding  Nj^  is  in 
a  direction  to  cause  the  core  to  set.  If  the  core  were  initially  set,  the  applica¬ 
tion  of  I2  through  N2  would  cause  the  core  to  reset.  The  magnetic  logic  in 
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connection  with  this  paper  is  performed  with  a  two-phase  clock.  The  core  is 
initially  in  the  reset  state.  At  phase  one  the  core  is  set  or  not  set  in  accordance 
with  the  logical  requirements.  At  phase  two  time  the  core  is  reset.  Signals 
induced  by  flux  change  are  interpreted  as  logical  ones. 

SET 


Figure  10  -  Mirror  Symbol  Notation  for  Core  Winding 

THE  VOTING  FUNCTION 

The  voting  function  is  a  Boolean  function  which  states  that  if  p  or  more  out 

of  n  variables  are  present,  then  an  output  will  result.  For  example,  the  truth 

table  for  p  =  2  and  n  =  3  is  shown  in  figure  11.  This  function  may  be  mechanized 

with  a  single  core.  It  is  assumed  that  a  unit  current  pulse  is  available  to 

56 

represent  each  variable  plus  its  prime.  Karnaugh's  rule  for  generating  the 
voting  function  is  to  connect  (p  -  1)  of  the  primed  variables  in  a  direction  to 
reset  the  core  and  (n  +  1  -  p)  unprimed  variables  in  a  direction  to  set  the  core. 
One  way  to  wire  the  core  is  shown  in  figure  12.  The  m.m.f.  equation  is 

•  '2'*2  ♦  ‘3“3  -  ‘l'”l 
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assuming  that  the  magnetic  threshold  of  the  core  is  one  ampere-turn  and  noting 
that  all  quantities  in  the  m.m.f.  equation  are  integers  it  may  be  stated  that  if 

F,  s  1  then  r.  =  1 
i  1 

and  F.  s  0  then  f.  =  0 

1  1 

where  the  index  "i"  denotes  the  row  of  truth  table  and  F^  and  f^  the  corresponding 
values  of  the  m.m.f.  and  Boolean  function  respectively.  It  is  readily  seen  that 
if  =  Ng  =  Ng  =  1  and  if  the  values  of  the  variables  are  applied  to  (13)  in 
accordance  with  the  truth  table,  the  proper  values  of  F^  are  attained  to  satisfy 
the  voting  function  required. 
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Figure  1 1  -  Truth  Table  for  Voting  Function 
where  p  =  2  and  n  =  3 
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Figure  12  -  Core  and  Wiring  Configuration  for  Voting 
Function  where  p »  2  and  n  =  3 
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It  should  be  noted  that  a  primed  variable  can  be  replaced  by  a  unit  current 
bias  pulse,  always  present  at  phase  one  time,  and  by  the  unprimed  variable  in 
accordance  with  (14). 


Complementing  is  achieved  since  I^' 
Equation  (13)  may  be  written  as 


I 

n 


=  1  when  I 

n 


(14) 

0  and  when  I  '  =  0,  I  =  1. 
n  n 


F  =  -IN  +  IN.  +  I„N,  +  I-N, 
oo  ll  a2  33 

where  N  has  the  same  number  of  turns  as  N. . 
o  1 

An  extension  to  this  type  of  switching  function  mechanization  is  to  develop 
a  technique  for  synthesizing  any  arbitrary  Boolean  function  in  a  similar  manner. 
It  turns  out  that  it  is  possible  to  mechanize  many  Boolean  functions  with  a  single 
toroid  appropriately  wound.  Unfortunately,  not  all  functions  can  be  realized 
with  a  single  core.  For  example,  the  "exclusive-OR"  requires  two  cores.  How¬ 
ever,  the  method  to  be  described  will  handle  multicore  solutions. 

Consider  the  function  01011  10^  (^denotes  least  significant  bit).  The 
truth  table  representation  is  shown  in  figure  13 .  If  the  core  is  wound  according 
to  (15),  which  means  I  and  I 

O  2 

F  =  21  +  I„  -  21.  -  I-  (15) 

o  z  13 

applied  to  two  turns  and  one  turn  respectively  in  a  direction  to  set  the  core  and 
I.^  and  Ig  applied  to  two  turns  and  one  turn  respectively  in  a  direction  to  reset 
core,  then  the  function  is  mechanized  in  a  single  core.  This  solution  may  be 
reached  by  intuition  or  by  trial  and  error.  When  the  number  of  variables 
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increases,  the  need  for  an  algorithm  becomes  apparent. 


THE  SIMPLEX  ALGORITHM 

A  linear  programming  technique  referred  to  as  Simplex  provides  an  algorithm 
for  determining  the  wiring  configuration  for  mechanizing  any  arbitrary  Boolean 
function  with  a  single  core  where  possible  and  with  a  multiplicity  of  cores  when 
necessary.  It  is  the  aim  of  this  paper  to  provide  an  insight  to  the  workings  of 
Simplex  and  its  application  to  our  problem,  while  at  the  same  time  avoiding  the 
many  associated  details  in  connection  with  computation  and  programming.  Con¬ 
sider  the  truth  table  of  figure  13.  An  m.m.f.  inequality  may  be  written  for  each 
row.  A  threshold  value  of  one  ampere-turn  is  assumed  for  core  switching,  and 
in  conditions  where  core  switching  is  not  desired  the  net  m.m.f.  is  held  less 
than  or  equal  to  zero  ampere-turns. 
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Figure  13  -  Truth  Table  for  f  »  0101 1101 
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loNo  +  ^  0 


I  N 

0  0 

'3N3 

(16) 

I  N  +  I^N, 
00  11 

^  0 

IN  + 

0  0 

I2N2  +  I3N3 

I  N  +  I,N,  + 
00  11 

I2N2  +  I3N3 

<,0 

The  system  of  inequalities  (16)  is  consistent  (one  or  more  solutions),  and 
the  function  may  be  mechanized  by  a  single  core.  If  the  system  of  inequalities 
(16)  were  inconsistent,  then  more  than  one  core  would  be  required. 

If  an  arbitrary  Boolean  function  is  processed  by  Simplex,  the  results  are 
as  follows: 

1.  If  the  system  of  inequalities  is  consistent,  a  set  of  values  for  will 
be  found  which  satisfy  the  inequalities.  However,  it  is  important  to  note  that 
a  particular  solution  will  be  found  which  will  minimize  an  arbitrarily  specified 
linear  function  referred  to  as  the  objective  or  cost  function  of  the  form 

Z  =  K,  N,  +  K„  N„  +  .  .  .  +  K  N  (17) 

1  1  «  6  n  ri 

where  K  are  constants.  (In  other  words  the  Simplex  algorithm  furnishes  a 

method  for  solving  a  system  of  inequalities  and  at  the  same  time  provides  a 

minimization  criterion,  namely  the  total  number  of  turns.  ) 

2.  If  the  system  of  inequalities  is  inconsistent,  then  the  Simplex  solution 
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will  indicate  how  the  constraints  must  be  modified  and  what  portions  of  the  desired 

function  will  be  mechanized  in  each  of  the  multiplicity  of  cores. 

Before  proceeding  with  the  explanation  of  the  Simplex  algorithm,  examples 

of  solutions  of  two  functions,  a  single  core  solution  and  two  core  solution,  are 

given.  The  solution  to  the  function  0101110^  is  =  2,  =  -  2,  N2  =»  1»  and 

N„  =  -  1.  This  of  course  results  in  (15).  Consider  the  function,  h  =  0110001£. 
o  *” 

The  Simplex-derived  solution  is  N,  =  1,  N„  =  -  1,  and  d  =  K.  Noting  the 
presence  of  d.  =  K  in  the  solution  and  temporarily  ignoring  it,  the  m.m.f. 

O 

equation  may  be  written. 

F  =  Ij  -  I2  (18) 

A  core  wound  in  accordance  with  (18)  yields  the  function  f  =  0010001^.  Note 
that  f  differs  from  h  by  its  absence  of  a  "one"  in  the  6th  position.  This  corres¬ 
ponds  to  the  subscript  of  d,  a  dummy  variable,  the  role  of  which  will  be  described 
later. 

A  new  function,  g,  is  formed  by  observing  the  dummy  variables  remaining 
in  the  solution  and  placing  a  "one"  in  the  bit  positions  corresponding  to  the  sub¬ 
scripts  of  the  dummy  variables.  In  this  example  g  =  0100000£.  This  function  is 
now  processed  by  Simplex,  and  the  solution  is  N  =  -  1,  N  =  -  1,  N„  =  1,  N_  *  1. 
The  corresponding  m.m.f.  equation  is 

G  =  -I  -  I,  +  I„  +  I_  (19) 

o  1  z  J 

Since  no  dummy  variables  remain  in  the  solution  to  g,  it  is  a  one-core  function. 

In  general  g  may  not  be  a  single-core  function,  and  the  procedure  described  above 
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would  be  repeated  until  the  final  function  processed  contained  no  dummy  variables. 
In  our  example,  two  cores  are  required  to  mechanize  h  (refer  to  figure  14).  Cores 
F  and  G  are  wound  according  to  (18)  and  (19)  respectively.  An  output  wire  is 
series  connected  to  both  cores.  Either  core  F  or  G,  or  neither,  will  switch  in 
accordance  with  the  input  lines  energized.  Hence,  the  voltage  induced  in  the 
output  line  at  phase  two  time  represents  the  function  h,  or 

(20). 


OUTPUT*  h 


Iz- 

u- 

RESET- 


-5t 


Figure  14  -  Two-Core  Mechanization 
of  h  =  011000010 
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EXPLANATION  OF  SIMPLEX 


Consider  the  following  system  oi  inequalities. 


(1) 

^1^1  +  ^2*2 

(2) 

^21*1  ^22*2 

(3) 

^31*1  ^32*2 

^bs 

(4) 

^41^1  +  ^42*2 

^>4 

(21) 


If  each  inequality  of  (21)  is  treated  as  an  equation  and  is  graphed,  the  set  of  lines 
shown  in  figure  15  would  result.  Only  the  area  included  by  the  convex  polygon  abed 
contains  values  of  and  X2  which  are  solutions  to  (21).  A  useful  theorem  enables 
one  readily  to  find  a  solution  which  will  minimize  an  associated  cost  function. 


-44- 


Theorem  I;  A  linear  function  defined  over  a  convex  polygon  takes  on  its 

maximum  and  minimum  value  at  a  corner  point  of  the  convex 
polygon. 

Hence,  one  simply  tests  the  value  of  the  cost  function  at  each  corner  of  the 
convex  polygon  to  determine  the  valued  of  x  and  x  which  satisfy  the  system  of 
inequalities  and  minimizes  the  cost  function.  As  the  number  of  variables  and 
equations  increases,  the  geometrical  picture  becomes  hopeless.  It  would  amount 
to  finding  corner  points  of  convex  hulls  in  hyperspace.  Fortunately,  the  Simplex 
algorithm  need  not  be  considered  in  the  light  of  the  above.  The  algorithm  is 
perhaps  best  explained  by  example  accompanied  by  pertinent  theorems  and 
definitions. 

Consider  the  following  system  of  inequalities  and  cost  function. 


(1) 

*1 

0 

(2) 

^2 

0 

(3) 

Xi  +  2X2 

$ 

8 

(4) 

*1  *2 

5 

(5) 

- 

5 

6 

Z  =  -  2x^  -  SXg 

The  convex  polygon  is  shown  in  figure  16.  The  value  of  the  coordinate  of  the 
corner  point  (14,  15)  minimizes  the  cost  function. 

The  Simplex  algorithm  is  utilized  to  find  what  is  referred  to  as  a  basic 
feasible  solution. 

Definition  I.  A  feasible  solution  is  a  solution  which  contains  nonnegative  x's. 
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Definition  11. 


A  basic  feasible  solution  is  a  solution  which  contains  at  most 


m  positive  x's  where  m  is  the  number  of  constraints. 

Since  the  solutions  are  restricted  to  basic  feasible  solutions,  the  first  two  in¬ 
equalities  of  (22)  need  not  be  written.  The  last  three  inequalities  of  (22)  are 
converted  to  equations  by  the  introduction  of  "slack"  variables  and  Xg. 


Figure  16  -  Corner  Points  of  the  Convex 
Polygon  of  Equations  (22) 
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8 


Z 


*1  +  2X2  +  Xg 


Xi+  X2  +  X4 


-  +  2X2  +  Xg 


2^1  -  3^2 


5 

6 


(23) 


Equations  (23)  are  now  said  to  be  in  "cannonical"  form  with  respect  to  variables 
Xg,  x^,  and  Xg.  Notice  that  the  coefficient  of  Xg  is  plus  one  and  that  it  does  not 
appear  in  any  other  equation.  This  is  true  also  for  x^  and  Xg.  The  canonical 
form  in  general  is 


^1,  m+  1  *m+  1  ^  ' 

.  .  +  a.  X 

In  n 

=  *^1 

^2,  m+  1  *m+  1  ^  ' 

.  .  +  a„  X 

2n  n 

II 

^n,  m+  1  *m+  1  ^  ' 

.  .  +  a  X 

mn  n 

*  b 

m 

C  ^  ,  X  +  .  . 

m+1  m+1 

.  .  +  C  X  + 

n  n 

constant 

(24) 


Theorem  II:  If  a  problem  has  an  optimum  feasible  solution,  then  it  also  has 

a  basic  feasible  solution  which  is  optimal. 

Note  that  the  canonical  form  of  (23)  enables  one  to  obtain  a  basic,  feasible 
solution  simply  by  setting  Xj^  and  equal  to  zero.  One  solution,  then,  is  Xg  =  8, 

=  5,  Xg  =  6,  and  the  corresponding  value  of  the  cost  function  is  Z  =  0.  The 
introduction  of  slack  variables  provides  an  initial  basic  feasible  solution  not 
necessarily  optimum.  By  appropriate  algebraic  manipulation,  (23)  can  be  put  into 
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Ceinonical  form  with  respect  to  any  set  of  three  x’s.  These  variables  are  referred 
to  as  "basic  variables"  or  as  the  "basis. "  In  this  problem  the  total  number  of 
solutions  to  be  examined  is  the  combination  of  five  things,  taken  three  at  a  time, 
or  ten  solutions.  For  larger  problems  this  trial  and  error  procedure  becomes 
prohibitive.  Fortunately,  there  is  a  systematic  iterative  procedure  for  finding, 
in  a  relatively  few  number  of  steps,  a  solution  which  minimized  the  cost  function. 

Referring  to  the  cost  function  of  (23),  if  x^^  or  x^  were  increased  from  zero, 
then  the  value  of  Z  would  be  reduced.  This,  of  course,  is  because  their  co¬ 
efficients,  called  the  "relative  cost  factors,  "  are  negative.  First  to  be  considered 
is  X2,  since  it  would  appear  that  any  change  here  would  have  the  greater  effect 
on  the  value  of  Z.  (However,  it  is  not  essential  that  x  be  given  first  considers- 

A 

tion.  )  The  question  arises  as  to  how  large  x  can  be  made  without  causing  any 

A 

of  the  other  x's  to  turn  negative.  Therefore,  from  (23)  the  following  may  be 
written: 


The  minimum  positive  value  of  is  selected;  i.  e.  Xg  =  3.  If  the  value  of 

4  were  selected,  then  x  would  turn  negative.  If  the  value  of  5  were  selected, 

0 

then  both  x  and  x  would  turn  negative.  The  value  of  x  =  3  makes  x  =  0. 

u  0  2  v 

Hence,  x  replaces  x  in  the  basis  and, the  equations  (23)  are  now  rearranged  to 
2  0 

canonical  form  with  respect  to  the  basic  variables  X2,  x^,  x^. 
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2Xi  -  Xg  +  Xg 


2 


(25) 

3 


The  basic,  feasible  solution  is  x^  =  0,  x^  =  3,  Xg  =2,  x^  =  2,  and  Xg  =  0 
with  the  cost  function  decreasing  to  Z  =  -  9.  By  examining  the  relative  cost 
factors  of  (25),  it  can  be  seen  that  the  value  of  the  cost  function  can  be  decreased 
by  increasing  x^^  only.  The  positive  value  of  the  relative  cost  factor  of  Xg  prohibits 
increasing  x^  from  zero  since  this  would  result  in  an  increase  for  Z.  Testing 
again  to  find  the  new  basis,  it  can  be  seen  that  x^  replaces  Xg. 


Notice  that  one  value  for  x^^  is  -6.  This  cannot  be  used  since  it  would  result 
in  an  increase  in  the  value  of  the  cost  function.  Equations  (25)  are  now  put  in 
canonical  form  with  respect  to  the  basic  variables  x^^,  Xg,  and  x^. 
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1 


^3  ^5 

2  “  2  ^1  “ 


3*3  Xg  1 

~  T'*'  ^4  “  T" 


(26) 


^3  ""S  7 

T~  +  *2  ““T 


Z  =  - 


^5  '^^3  25 

4  ■’■  4  “2 


The  solution  is  =  1,  =  2  *  *3  “  ^4  ~  2  »  *5  “  value  of 


„  .  25 

Z  IS  -  2  •  Repeating  the  process 


=  1 


X5  =  -2 


-i  .  i_ 
4  ■  2  ’ 


X5  =  2 


j__  7 
4  "  2  ’ 


X5  .  14 


The  new  basis  is  x,,  x„,  and  x_ 
i  J  5 


*3  ■*■  3x^  +  Xj  =  2 


3x3  +  4x^  +  Xg  -  2 


(27) 


X3  -  X4  +  x^  .  3 


z  *  -^Xg  +  x^  -  13 
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The  solution  is  x„  =  0,  x.  =  0,  x_  =  2,  x  *  2,  and  x„  =  3.  The  value  that  Z  takes 
j  4  a  1  z 

is  -  13.  Note  that  now  all  of  the  relative  cost  factors  are  positive.  Hence,  neither 
Xg  nor  x^  can  be  increased.  This  indicates  that  the  above  solution  is  one  which 
minimizes  Z.  This  solution  agrees  with  the  solution  obtained  by  the  graphical 
method.  Note  that  the  intermediate  solutions  are  also  corners  of  the  convex 
polygon  of  figure  16. 

Since  the  Simplex  solution  yields  only  positive  values  for  the  unknowns,  each 

N.  is  represented  as  the  difference  of  two  positive  numbers  (t.  -  C.)  which  would 
1  J  3 

correspond  to  clockwise  and  counter-clockwise  windings,  respectively.  Hence, 
the  solution  to  function  01011101  would  be  of  the  form  t  =  2,  C  =  2,  t„  =  1,  and 
Cg  =  1.  Preparing  inequalities  (16)  for  Simplex  they  would  be  written  as 


o  o 


O  0 

+ 

2 

0 

0  o 

+ 

1 

o  o 

+ 

+ 

1 

0  o 

+ 

1 

(28) 

) 

o  o 

+ 

(ti-Ci) 

+ 

0 

(t -c^) 
o  0 

+ 

+ 

<^3-S) 

1 

(t  -C  ) 
o  o 

+ 

+ 

<^2-^2) 

+ 

2 

0 
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where  the  t.  and  C.  are  the  unknowns.  The  I  coefficients  are  dropped  because 
J  J  n 

they  have  the  value  of  unity. 

Note  that  in  addition  to  s  signs  appearing  in  (28)  ^  also  appear.  In  the 
case  of  ^  constraints  a  slack  variable  is  subtracted  and  a  dummy  variable  is 
added  to  the  lefthand  side  of  the  inequality.  Dummy  variables  are  required  to 
provide  the  initial  basis.  The  initial  basis  is  then  comprised  of  slack  Eind  dummy 
variables,  namely  S^,  d^,  d^,  d^,  S^.  S^,  d^,  d.^. 


o  o 


+  S 


(t  -  C  )  +  (t  -  CJ 
o  o  1  1 


K  - 

o  o 


^  <^2  -  ''2> 


(t  -  C  )  +  (t,  -  C.)  +  (t„  -  C„) 

O  O  11  11 


d^=  0 

^2  ^  ‘^2  =  ' 
S3  ^3  ^  ^ 


(29) 


<^o  - 

<*o  -  ^o)  "  "  ^1> 


+  (t3-C3)+S5  =0 


o  o 


+  (t2  -  c^)  +  (tg  -  C3)  -  Sg  +  dg  =  1 


<*o  -  ''o)  +  <*1  -  C,)  +  (t2  -  C^)  4.  (tg  -  C3)  -  4  d^  =  0 


The  cost  function  is 

Z  =  t  4  t,  4  t„  4  t„  4  C  4  C,  4  C„  4  C„  4  Md,  4  Md„  4  Md.  4  Md„  4  Md„  (30) 

01«001aiJ  1  1  O  O  I 

Note  that  dummy  variables  are  included  in  the  cost  function.  The  coefficients  M 
of  the  dummy  variables  are  chosen  very  large  to  bias  the  problem  so  that,  if  possible, 
the  dummy  variables  will  not  be  included  in  the  basis  which  minimizes  Z*, 

In  some  Simplex  applications  the  existence  of  a  dummy  variable  in  the  final  basic 
indicates  there  is  no  feasible  solution.  In  our  application  this  type  of  result  sim  ply 
indicates  that  more  than  one  core  is  required  to  mechanize  the  given  function. 
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Equation  (30)  is  rearranged  to  canonical  form  by  appropriate  substitutions 


from  (29),  eliminating  the  dummy  variables. 

Z  =  (1  -  5M)t^  +  (1  -  3M)t^  +  (1  -  4M)t2  +  (1  -  2M)tg  +  (1  +  5M)C^ 

+  (1  +  3M)C^  +  (1  +  4M)C2  +  (1  +  2M)C3  +  MS^  +  MS2  +  MSg  +  MSg  +  MS^ 

+  5M 

The  problem  is  now  ready  to  be  treated  in  the  manner  previously  described. 

It  will  be  recalled  that  if  a  function,  h,  is  processed  by  Simplex,  and  dummy 
variables  appeared  in  the  solution,  then  a  multi-core  configuration  was  required 
to  mechanize  the  function.  Further,  if  one  core  were  wound  in  accordance  with 
the  F  and  resulting,  part  of  the  function  was  obtained  and  a  new  function  was 
defined  by  placing  a  "one"  in  positions  indicated  by  the  subscripts  of  the  dummy 
variables  in  the  solution.  Let  us  refer  to  the  portion  of  the  function  initially 
mechanized  as  f,  and  the  remaining  portion  as  the  residue  function,  g.  Equation 

(29)  may  be  written  as  follows: 

when  h^  =  1  then  ^  ^  (t^  -  C^)  -  S.  +  d.  =  1 

h.  =  0  then  V]  .  (t  -  C  )+  S.  *  0 

Let  F.  =  Z.  (t.  -  C.);  hence,  when  F.  i  1,  f.  =  1,  and  F.  s  0,  L  =  0. 

i  3  J  J  ^  ^ 

The  symbol  [  ]  defined  as  follows: 

[x]  =1  if  X  i  1 

[x]  .0  if  X  ^  0 


(31) 
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The  residue  function  is  defined  as  follows: 


1 

0 


when 

when 


d.  >  0 
1 

d.  =  0  or  when  there  is  no  d. 

1  1 


Equation  (20)  may  now  be  written  as 

^  =  [•'i]  *  [°i] 


To  prove  (34)  all  possible  combinations  will  be  examined;  i.  e,  : 


h.  = 
1 


h. 

1 


=  1 

=  0 

=  1 

=  0 


1,  h.  .  l;  [fJ  .  1 


i.e.  F.  i  1 
1 


(32) 

(33) 


(34) 


(35) 


Since  Simplex  guarantees  that  the  original  set  of  equations  must  be  satisfied,  we 


may  write 

F.  -  S.  +  d.  =  1 
111 

But  from  (35) 

F.  i  1 

1 

(36) 

also 

Sj  i  0 

(all  variables  are  positive) 
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Adding  (36) 


2F.  +  d.  i  2 
1  1 


or 


(37) 


Hence  where  ®  1  i  1  the  dummy  variable  must  take  the  value  of  zero 

(all  variables  are  positive)  and  we  may  say  from  (33)  that  =  0. 


2.  h.  =  1; 


=  0 


1.  e. 


F. 


F.  -  S.  +  d.  =1 
1  11 


(38) 


-F. 

1 


S. 


^  0 
^  0 


Adding  (38) 


d.  ^  1  or  d.  >  0 

1  1 


=  0 

1 


Hence,  where  /  F.  |  =  0  dummy  variables  must  remain.  From  (32) 


w- 


3,  h.  =  0; 
1 


[•'i] 


=  1 


i.e.  F.  ^  1 
•  1 


F.  +  S.  =0 
1  1 


(39) 


In  this  case,  (39)  could  be  satisfied  only  if  were  negative,  which  is  not  permitted 


by  Simplex.  Hence,  the  condition  (15)  is  not  possible. 


4.  h.  =  0;  fpr]  =0  i.e.  F.  s  0 

^  Lil 

In  this  case  (39)  is  satisfied  by  a  positive  slack  variable  which,  of  course,  is 
permissible.  However,  since  there  is  no  d^,  we  may  say  from  (33), 
Therefore,  equation  (34)  is  satisfied  under  all  possible  combinations. 


0. 
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The  residue  function  defined  above  is  now  treated  by  Simplex.  If  it  is  a 
one-core  function,  the  problem  is  terminated,  yielding  a  two-core  function.  If, 
in  processing  the  residue  function,  dummy  variables  remain  in  that  solution, 
then  a  new  residue  function  is  defined. 

The  above  procedure  is  repeated  until  the  final  residue  function  processed 
contains  no  dummy  variables.  Therefore,  one  magnetic  core  is  required  for 
each  function  or  residue  function  processed  in  the  course  of  a  problem.  From 
another  point  of  view  it  may  be  said  that  the  given  function  in  a  problem  is  re¬ 
duced  to  a  set  of  one -core  functions,  each  of  which  might  be  represented  by  a 
consistent  set  of  inequalities.  To  obtain  the  original  function,  all  of  these  one- 
core  functions  are  ORed  together  by  the  output  line. 

SIMPLEX  MODIFICATIONS 

In  the  course  of  the  applications  of  the  Simplex  algorithm  to  the  above 
problem,  several  difficulties  arose  which  resulted  in  some  modifications.  For 
example,  some  solutions  yielded  minimum  total  number  of  turns  at  the  expense 
of  additional  cores.  This,  under  the  requirements  of  unlimited  total  number  of 
input  turns,  was  not  desirable.  This  type  of  answer,  however,  may  be  quite 
useful  and  will  be  discussed  under  future  work.  An  extreme  example  of  this  type 
of  difficulty  is  the  result  calling  for  zero  turns  and  a  residue  function  equal  to 
the  original  function.  The  type  of  Simplex  algorithm  finally  developed  is  "dummy- 
oriented;"  i.  e.,  m  every  applicable  decision  process  during  the  course  of  the 
problem,  dummy  variables  are  removed  from  the  basis.  The  following  modifi¬ 
cations  were  introduced; 
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1.  Once  a  dummy  leaves  the  basis  it  is  not  permitted  to  return.  This  is 
accomplished  by  ignoring  the  relative  cost  factors  of  dummy  variables  in  the 
cost  function. 

2.  The  problem  is  not  necessarily  permitted  to  end  when  all  of  the  relative 
cost  factors  become  positive.  Instead,  the  initial  conditions  of  the  problem  are 
altered  by  changing  the  value  of  a  relative  cost  factor  assigned  to  one  of  the 
dummy  variables  so  that,  at  this  stage  of  the  problem,  one  of  the  non-dummy 
variables  assumes  a  negative  relative  cost  factor.  The  test  described  previously 
is  applied  to  determine  which  variable,  if  any,  will  leave  the  basis.  It  is 
possible  that  a  variable  with  a  negative  relative  cost  factor  cannot  enter  the  basis 
because  the  test  for  minimum  positive  value  yields  only  negative  numbers.  It 
should  be  noted  that  it  is  not  difficult  to  change  the  initial  conditions  in  the  manner 
described  above  because  only  the  cost  equation  is  affected  throughout  the  problem. 

This  procedure  may  be  explained  from  the  geometrical  point  of  view.  Con¬ 
sider  the  cost  function  of  (22).  If  the  cost  function  were  changed  to  Z  *  -  SOx^, 

then  a  different  corner  of  the  convex  polygon  would  minimize  the  new  cost  function. 

Summarizing,  it  may  be  stated  that  the  relative  cost  factors  are  repeatedly 
altered,  forcing  the  problem  to  continue  until  a  basis  containing  a  minimum  number 
dummy  variables  is  obtained. 

At  the  price  of  a  more  "expensive"  solution,  i.  e.,  one  which  might  contain 
more  turns  but  uses  less  cores,  the  problem  is  permitted  to  terminate  only  when 
it  is  impossible  to  remove  dummy  variables  from  the  basis. 

It  is  interesting  to  note  that  it  is  always  possible  to  remove  at  least  one 
dummy  from  the  initial  basis.  The  partial  result  would  be  a  function  containing 
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a  single  "one"  and  a  residue  function  with  one  less  bit  than  the  original  function. 
All  Boolean  functions  containing  a  single  "one"  are  realizable  with  a  single  core. 
But  once  a  particular  dummy  leaves  the  basis,  all  dummy  variables  in  rows 
corresponding  to  ones  which  were  consistent  with  the  first  dummy  to  leave  the 
basis  will  also  leave  the  basis. 

VARIABLE  SUPPRESSION 

A  method  referred  to  as  "variable  suppression"  has  been  developed.  This 
method  enables  the  synthesis  of  n-variable  functions  in  terms  of  a  known  solution 
of  functions  of  less  than  n  variables.  The  variable-suppression  technique  yields 
solutions  which  may  not  be  minimal,  but  this  technique  is  one  which  is  readily 
automatized. 

This  technique  is  best  described  by  example.  For  simplicity,  a  function 

of  four  variables  will  be  synthesized  in  terms  of  known  solutions  of  functions  of 

two  variables.  Consider  the  function  f  =  110100101  lOOlOH.  The  truth  table  is 

shown  in  figure  17.  Note  that  the  function  f  is  divided  into  four  subfunctions,  a,  b, 

c,  and  d.  The  portions  of  each  function  lined  out  by  arrows  are  zeros.  It  will 

be  observed  that  f.  =  a.  +  b.  +  c.  +  d..  Hence,  the  synthesis  of  the  subfunction 
11111 

will  provide  a  solution  for  f  simply  by  connecting  the  output  winding  of  the  cores 

associated  with  a,  b,  c,  and  d  in  series. 

Consider  function  a.  Bit  positions  i  =  0  through  i  =  3  may  be  considered  as 

a  function  of  two  variables,  namely  I  and  I  The  solution  for  this  function  is 

X  « 

I  +  I,  -  I„  (36) 

O  \  Z 
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1,1. 1. 1. 

f 

1  t 

( 

1*0 

0  0  0  0 

'  ' 

1 

0  0  0  1 

• 

2 

0  0  10 

0  0 

S 

0  0  11 

1 

4 

0  10  0 

0 

> 

S 

0  i  0  ( 

0 

0 

6 

0  110 

1 

7 

0  111 

1 

6 

10  0  0 

0 

0 

9 

I  0  0  i 

1 

1  0 

10  10 

0 

0 

1  1 

(Oil 

0 

0 

1  2 

110  0 

1 

1 

1  S 

1  1  0  i 

0 

0 

1  4 

1  1  t  0 

1 

1 

1  9 

1  (  i  ( 

1 

1 

Figure  17  -  Truth  Table  for  f  and  the  Sub  Functions 

But  it  is  not  sufficient  to  wind  core  "a"  with  this  input  configuration  because  it 
represents  a  function  of  two  variables  which  would  repeat  itself  every  four  bit 
positions.  It  is  necessary  to  suppress  any  ones  that  might  appear  in  the  re¬ 
maining  portion  of  the  function. 

Note  that  the  most  positive  value  (36)  can  take  is  2.  If  a  winding  of  two  turns 
is  provided  for  variable  I  and  negatively  connected,  then  the  presence  of  the  I 
current  pulse,  during  bit  positions  i  =  4  through  i  =  7,  would  cause  zeros  in  these 
positions.  Similarly,  a  winding  for  of  two  turns  negatively  connected  would 
cause  zeros  in  bit  positions  i  =  8  through  i  »  15.  Hence,  the  m.m.f.  equation 
for  "a"  is 

A  =  X  +  X  -  X.,  -  2X,  -  2X, 
o  1  2  3  4 
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where 


a.  =s  1 

1 

when 

A.  i 

I 

1 

a.  =  0 

1 

when 

A.  ^ 

1 

0 

The  solution  for  bit  positions  i  =  4  through  i  =  7  for  function  b  is  simply  I  . 

Bit  positions  i  =  0  through  i  =  3  are  held  to  zero  by  one  turn  negatively  connected 
for  I  The  equivalent  of  I  '  is  I  -  I~.  Bit  positions  i  =  8  through  i  ■  15  are 

O  o  O  o 

held  to  zero  by  one  turn  negatively  connected  for  variable  X^.  The  m.m.f. 
equation  for  b  is 

B  =  X  +  X-  +  X„  -  X . 
o  2  3  4 

Functions  c  and  d  are  handled  in  the  manner  described  above.  Summarizing,  the 
set  of  m.m.f.  equations  for  mechanizing  f  are 

=  ^o  ^1  -  ^2  ‘  2X3  -  2X4 
B  =  -X  +  X„  +  X„  -  X . 

0^34 

C  =  -X^  +  X  -  X,  -  X„  +  X. 
o  12  3  4 

D  =  -3X  -  X,  +  X.,  +  2X  ‘  +  2X^ 
o  i  2  3  4 

The  cores  and  the  winding  arrangement  are  shown  in  figure  18. 

The  Simplex  algorithm  previously  described  was  programmed  and  operated 
successfully  on  the  Burroughs  220  Computer.  This  technique  which  will  be 
referred  to  as  Simplex  Mod  I  provided  the  capability  of  handling  switching  functions 
of  six  variables.  Computer  storage  capacity  is  the  limiting  factor  with  respect 
to  the  number  of  variables  which  may  be  processed. 
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OUTPUT  =  f 


Figure  18  -  Cores  and  Winding  Arrangement  for  Mechanizing  f 


An  improved  version  of  Simplex  (Simplex  Mod  II)  has  also  been  programmed 
for  the  220  which,  in  addition  to  reducing  computation  time,  has  the  capability 
of  handling  switching  functions  of  nine  variables.  The  modifications  made  to  the 
former  method  tffect  the  tableau  structure,  the  computational  procedure,  and 
the  cost  function,  and  are  described  below. 
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TABLEAU  STRUCTURE 


It  will  be  recalled  that  the  constraining  equations  are  the  m.m.f.  inequalities 

written  for  each  row  of  the  truth  table.  In  order  to  provide  for  both  clockwise 

and  counterclockwise  turns,  each  is  represented  as  the  difference  of  two 

positives  (t.  -  c.)  which  correspond  to  clockwise  and  counterclockwise  turns 
J  J 

respectively.  Consider  the  two-variable  function  f  =  1001.  The  inequalities  are: 


(^o  - 

K  ^ 


^  1 
s  0 
^  0 


(t  -  c  )  +  (t  -  c  )  +  (t  -  c  )  i  1 
00  11  2  2 


where  t.  and  c.  are  the  unknowns. 


(37) 


J  3 

In  forming  an  equation  from  an  inequality  containing  the  ^  sign,  a  slack 
variable  is  added.  For  those  inequalities  containing  the  2  sign  a  slack  variable 
is  subtracted.  At  this  point  there  is  a  deviation  from  Mod  I  in  that  a  dummy 
variable  is  not  also  added  to  the  left  side  of  the  equation.  Setting  up  the  constraints 
in  this  manner  permits  the  of  the  equations  with  2  signs  to  take  on  negative 
values.  Now  the  initial  solution  is  a  basic  solution,  rather  than  a  basic  feasible 
solution.  Variables  with  negative  values  are  permitted  and  given  physical 
interpretation. 

Considering  this  it  follows  that  the  initial  basis  is  comprised  of  slack  variables 

S  ,  S,,  S„,  S-  with  constraining  equations  of  the  following  form: 
o  1  2  3 
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1 


K  - 

O  O 

-  S 

( 

(t  -  c  )  + 
o  o 

(ti  -  C^) 

+  s 

(t  -  C  ) 

0  o 

+ 

<*2  -  ‘^2) 

+  s, 

+ 

0  o 

(t^  -  c^)  + 

^*2  -  ^^2) 

-  s, 

(38) 


In  the  Mod  I  a  cost  function  is  specified.  In  Mod  II  a  cost  function  is  selected 
from  one  of  the  constraints  which  is  of  the  form 


S.  =  EN.  -  1  (38a) 

1  J 

Where  the  summation  represents  the  total  number  of  input  turns  (disregarding 
polarity)  in  the  magnetic  core  and  the  minus  one  is  the  initial  value  of  the  slack 
variable,  the  quantity  to  be  maximized.  Only  those  slack  variables  with  a  value 
of  minus  one  are  eligible  to  be  selected  as  cost  functions.  Constraints  which 
are  in  contradiction  will  take  on  a  value  of  s:  -  2. 

Consider  the  example  f  *  1001.  The  inequalities  are 


X 

0 

1 

(39) 

X  +  X, 

0  1 

0 

(40) 

^0  ^2 

0 

(41) 

Xo  +  +.  X2 

1 

(42) 

adding  (39)  and  (42),  (40)  and  (41),  we  get 
2X^  +  i  2 

contradiction 

2X  +  X,  +  :£  0 

o  1  2 
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Adding  the  necessary  slack  variables  to  the  inequalities  we  have: 


X 

0 

-  s  = 

0 

1 

(43) 

X 

0 

^  ^1 

0 

(44) 

X 

0 

+  +  Xg 

^  ^2  = 

0 

(45) 

X 

0 

+  X^  +  Xg 

-  S3  = 

1 

(46) 

Adding  (43)  and  (46),  (44)  and  (45-),  we  get 

2X  +  X,  +  X-  -  S  -  S-  =  2  (47) 

O  1  A  Q  i 

2X^  +  X^  +  X2  +  +  S2  =  0  (48) 

subtracting  42  from  41 

But  S2  and  are  initially  positive  quantities  and  are  not  permitted  by  the  Simplex 
Method  to  take  on  negative  values.  Removing  S2  and  Sg  from  (49)  we  get 

S  +  S,  s  -2  (50) 

O  O 

If  either  constraints  (43)  or  (46)  is  optimized  making  S  or  S„  ^  0,  the 

o  o 

burden  of  satisfying  (50)  is  left  to  the  other  giving  in  each  case  ^  -2  or 

Sg  ^  -2,  verifying  that  the  form  of  the  cost  function  must  be  as  previously  stated. 

At  first  glance  it  might  appear  that  in  maximizing  the  value  of  we  are  being 
forced  to  maximize  since  this  is  the  only  variable  term  in  the  cost  function. 

It  is  obvious  that  this  is  undesirable,  and  it  can  be  shown  that  not  only  can  be 
maximized  without  also  maximizing  but  that  the  minimum  value  of  EN^  is 

achieved.  Before  this  is  shown,  however,  it  might  be  best  to  ex  lain  why 
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is  being  maximized  and  why  the  objective  function  takes  the  form  of  (38a). 

It  will  be  recalled  that  if  a  function,  h,  is  processed  by  Simplex,  and  dummy 
variables  appeared  in  the  solution,  then  a  multi-core  configuration  was  required  to 
mechanize  the  function.  Further,  if  one  core  were  wound  in  accordance  with  the 
t^  and  resulting,  part  of  the  function  was  obtained  and  a  new  function  was 
defined  by  placing  a  "one"  in  positions  indicated  by  the  subscripts  of  the  dummy 
variables  in  the  solution.  In  Mod  II  dummies  are  replaced  by  negative  slack 
variables.  Refer  to  the  portion  of  the  function  initially  mechanized  as  f,  and 
the  remaining  portion  as  the  residue  function,  g.  The  constraints  may  be 
written  as  follows: 

when  h.  =  1  then  E(t.-c.)  -  S,  =  1  (51) 


h.  =  1  then  E(t.  -  c.)  -  S,  =  1 
1  J  J  1 


h.  =  1  then  E.(t.-c.)+S,=  0 
^  3  J  J  i 


(52) 


LetF  =  E.(t  -  c  )  ;  hence  when  F.  i  1,  f,  *  1  and  when  F.  s  0,  f.  =  0. 
1JJ3  11  11 

The  symbol  [  ]  is  defined  as  follows: 

[x]  =1  if  X  i  1 
l^x]  =0  if  X  s:  0 

The  residue  function  is  defined  as  follows: 

J  =  1  when  <  0 

^  G .  J  =  0  when  ^  0  or  when  there  is  no 

The  function  h  is  formed  by 

h.  =  [Fi]  +  [Gi] 


(53) 

(54) 

(55) 
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To  prove  (55)  all  possible  combinations  will  be  examined,  i.  e. , 

h.  .  1  [f,] 

=■  1 

h] 

■  0 

hi  =  0  [Fi] 

I  =  1 

hi 

=  0 

1.  h.  .  1,  [f.] 

=  1  i.e. ,  F.  s 

1  (56) 

Since  Simplex  guarantees  that  the  original  set  of  equations  must  be  satisfied. 

we  may  write 

F.  -  S.  =  1 

1  1 

(57) 

But  from  (56)  ”  ^  -1 

Adding  (57)  -  S.  s  0 

(58) 

Hence,  where  h.  *  1  and  F.  a  1,  the  slack  variable  must  be  i 

1  1 

0  and  from  (54), 

[Gi]  ■  0. 

2.  h,  =  [fi]  =  0 

i.  e. ,  .  F.  s  0 

1 

F.  -  S.  =  1 

1  1 

(59) 

o 

1 

adding  (59)  -S^  i  1  or  S^  <  0 

hence,  where  h^  =  1  and  j_FiJ  = 

0,  negative  slack  variables  must  remain. 

From  (53),  |_Gi]  =  I* 

3.  hi  .  0;  [Fi]  =  1 

i.e. ,  ^  1 

F.  +  S.  =  0 

1  1 
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(60) 

In  this  case.  (60)  could  be  satisfied  only  if  were  negative  which  is  not  permitted 
by  Simplex.  Hence,  the  condition  (39)  is  not  possible. 

4.  h^  =  0;  [F^]  =0  i.  e..  s  0 

F.  +  S.  =  0  (61) 

1  1 

-F.  ^  0 
1 

adding  (61) 

S.  ^  0  . 

1 

This  case  is  satisfied  by  a  positive  slack  variable  which  of  course  is  permissible. 
Since  S^  is  positive  we  may  say  from  (54),  [G^]  =  0.  Therefore,  equation  (55) 

is  satisfied  under  all  possible  combinations. 

The  residue  function  defined  above  is  now  treated  by  Simplex.  If  it  is  a  one- 
core  function,  the  problem  is  terminated,  yielding  a  two-core  function.  If,  in 
processing  the  residue  function,  negative  slack  variables  remain  in  the  solution, 
then  a  new  residue  function  is  defined. 

The  above  procedure  is  repeated  until  the  final  residue  function  processed 
contains  no  negative  slack  variables.  Therefore,  one  magnetic  core  is  required 
for  each  function  or  residue  function  processed  in  the  course  of  the  problem.  From 
another  point  of  view,  it  may  be  said  that  the  given  function  h  is  reduced  to  a  set 
of  one-core  functions,  each  of  which  can  be  represented  by  a  consistent  set  of 
inequalities.  To  obtain  the  original  function,  all  of  these  one-core  functions  are 
ORed  together  by  the  output  line. 

From  the  above  discussion  it  becomes  obvious  that  it  is  desirable  to  satisfy 

Slack  variables  whose  values  are  initially  positive  are  not  permitted  to  turn 
negative . 
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as  much  of  the  original  function  h  as  possible  by  successively  maximizing  as 
many  cost  fimctlons  (negative  slack  variables)  as  possible.  Looking  once  again 
at  the  cost  function 

S.  =  EN.  -  1 
1  1 

it  becomes  apparent  that  if  in  the  maximizing  process  is  made  zero,  the  corres¬ 
ponding  h^  is  satisfied.  Specifically  solving  the  first  equation  of  (38)  for  we 
get 

S,  =  (t  -  cj  +  (t  -  c  )  -1 
loo  11 

which  is  the  form  of  the  objective  function.  It  will  later  be  shown  that  the  maximum 

value  that  S.  can  take  is  zero. 

1 

Setting  up  equations  (38)  for  computational  procedure,  the  following  tableau 
is  produced: 


h 

i 

'  N 

0 

|!i 

N2 

S  ! 

o 

l 

|-1 

-1 

1 

To 

1 

0 

_  Si_  __ 

:  1 

1 

•  1 

0 

"2 

0 

1 

'0 

1 

\ 

:-l 

e 

-1 

-1 

For  purposes  of  tableau  interpretation  it  is  convenient  to  imagine  the  equations 
re-written,  so  that  all  variables  are  on  the  left  side  and  only  the  constants  on  the 
right  hand  side.  The  rows  refer  to  the  basic  variables  and  the  columns  to  the 
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non-basic  variables.  It  is  to  be  noticed  that  although  care  was  taken  in  setting 
up  the  constraint  equations  to  allow  for  positive  and  negative  solutions  by  repre¬ 


senting  the  N.  by  the  difference  of  two  positive  numbers  (t.  -  c.),  this  need  not 
J  J  J 

be  carried  into  the  tableau,  but  advantage  can  be  taken  of  the  fact  that  the  co¬ 
efficients  of  the  t.,  c,  pairs  in  the  non-basis  are  always  numerically  equal — but 
J  J 

opposite  in  sign.  Define  N.  =  t.  ;  N.  *  c..  This  definition  plus  the  removal  of  the 

J  J  J  J 

dummy  variables  of  Mod  I  results  in  a  dramatic  reduction  in  the  tableau  size 

for  Mod  II.  The  tableau  size  for  Mod  I  is  given  by  Rows  =  2^  and  Col  =  2(M=1)  +  K. 

M 

For  Mod  II  given  by  Rows  =  2  and  Col  =  M  +  1,  where -M  is  number  of  variables 
and  K  is  number  of  ones  in  the  function.  For  example,  a  nine-variable  problem 
containing  200  ones  in  the  function  requires  112,640  elements  in  the  tableau 
for  Mod  I  and  5, 120  elements  in  the  tableau  for  Mod  11. 


COMPUTATIONAL  PROCEDURE 

To  simplify  the  explanation,  and  to  arrive  at  a  clear  understanding  of  the 
modified  Simplex  process,  a  representation  sample  problem  is  examined  to 
demonstrate  the  computational  procedure.  For  reasons  of  computational  ease 
and  with  no  loss  of  generality  a  two-variable  function  with  a  multi-core  solution 
will  be  used;  specifically,  let  us  consider  h  s  1001.  The  constraining  equations 
are; 


N 

o 

-s 

0 

=  1 

* 

N 

o 

+ 

* 

^1 

=  0 

* 

N 

o 

+ 

♦ 

^2 

^  ^2  = 

♦ 

N 

o 

+ 

♦ 

+  N2 

-^3  * 

N*  =  (N.  -  N^) 

3  3  3 
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From  here  on  only  one  of  the  pair  will  be  carried. 

By  inspection  it  is  evident  that  a  basic  solution  is  given  by 

=  0.  =  0.  =  0.  =  -1.  =  0,  S2  =  0,  S3  =  -1  . 

The  constraints  for  the  basic  variables  are  expressed  in  terms  of  the  non-basic 
variables  and  the  initial  tableau  is  set  up  as  follows: 


cost 

,  ,  function 

1  ^  ' 

N 

0 

Ni 

^2 

(a)  S  =  N  -1  /  0 

0  0  L_w  S 

^  0 

-1  '-1 

0 

0  pivot  row 

(b)  Sj  =-N^  -Nj 

(^2)  -■ 

0  1 

1 

:  0 

S2=-N^  -N3 

A 

0  ;  1 

0  1 

(d)  S3  =  +  Nj  t  N3  -1 

0 

■-1  '-1 

-1  -1 

A  '• 

i  f pivot  col. 


Thus  far,  no  objective  function  has  been  determined,  and  before  proceeding  any 
further  one  must  be  picked.  Equation  (62a)  is  chosen  as  the  cost  function  because 
it  is  of  the  form  previously  described.  It  can  be  noted  that  equation  (62d)  is  also 
of  ttie  proper  form,  but  the  first  one  encountered  in  scanning  from  top  to  bottom 
is  always  picked. 

In  order  that  a  row  be  picked  as  the  cost  row  in  the  tableau  representation 
of  the  constraints  it  must  contain; 

(1)  A  slack  variable  in  the  cost  column  of  that  row. 

(2)  A  minus  one  in  the  corresponding  constant  column  element. 

(3)  A  negative  element  which  has  a  non-slack  variable  row  element. 
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At  each  step  one  of  the  non-basic  variables  is  brought  into  the  basis,  for 
which  the  entry  in  the  cost  function  has  a  negative  sign.  In  this  specific  problem 
at  this  point  is  to  be  made  basic.  The  column  headed  by  will  be  the  pivot 
column. 

It  should  be  noted  here  that  it  is  possible  and  quite  probable  that  there  will 
exist  more  than  one  non-basic  variable  which  would  have  these  conditions.  If 
row  (d)  is  picked  as  the  cost  row,  both  N  and  N„  would  be  eligible.  When  this 
situation  arises,  the  first  variable  encountered,  going  from  left  to  right,  is 
chosen. 

That  variable,  (or  one  of  those  in  case  of  a  tie)  which  has  the  smallest 
positive  ratio  of  the  constant  and  the  coefficient  in  the  column  thus  chosen,  must 
be  removed  from  the  basis;  only  slack  variables  are  taken  into  consideration 
when  determining  the  minimum  ratio.  This  row  is  the  pivot  row.  Again  looking 
at  the  tableau  will  leave  the  basis  since  row  two  has  a  ratio  of  zero. 

Going  back  to  the  constraining  equations  let  us  justify  all  that  has  been  said 
thus  far  about  the  choice  of  variables  for  maximizing  (62a).  can  only  be 
maximized  by  increasing  N^.  If  is  increased,  the  values  of  Sj^,  Sg,  and  Sg 
will  change.  Care  must  be  taken  at  this  point  not  to  make  either  or  Sg  negative. 
It  must  be  understood  that  in  this  problem  the  slack  variables  which  were  added 
to  the  inequalities  are  sign-restricted.  Therefore  they  must  remain  zero  or 
assume  positive  values.  Taking  this  into  consideration,  the  largest  value 
may  assume  is  zero.  In  the  attempt  to  maximize  S^,  no  change  is  effected  at 
this  time.  This  is  due  to  degeneracy  which  is  inherent  in  our  problem. 
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Since  N  was  selected  to  replace  S,  in  the  basis,  the  new  basis  will  consist 
o  1 

of  S^,  N^,  Sg,  Sg.  To  begin  the  next  phase,  these  variables  are  solved  in  terms 
of  the  non-basic  variables  getting  the  new  equations  below.  Using  the  rules  for 
tableau  transformation,  the  new  tableau  is  derived. 


h 

s 

1 

s 

s 

(a) 

S  =  -N, 
o  1 

- 

s 

0 

-1 

1 

1 

1 

-  ■  1 

0 

(b) 

- 

s 

(63) 

N 

0 

0 

i 

1 

1  i 

0 

(c) 

S2  =  +  Ni 

-Ng  + 

s 

s 

0 

-1 

-1 

1 

(d) 

S' 

+  Ng  - 

s 

-1 

1 

1  ® 

-1 

Equation  (63a)  no  longer  appears  in  a  form  suitable  to  pick  as  a  cost  function. 
Correspondingly;  the  first  row  of  the  tableau  fails  to  meet  the  requirements  of  a 
cost  row  because  it  does  not  contain  a  negative  element  with  a  non-slack  correspond¬ 
ing  variable  row  element.  It  must  be  remembered,  however,  that  in  order  to 

♦ 

reduce  the  size  of  the  tableau  only  one  part  of  was  chosen  to  appear  and 

recalling  also  that  N.  =  "N,.  Substituting  -N  for  N  ,  the  following  is  true 
J  J  1.1 


Cost 

Function 

h 

!■ 

s 

(a)  N,  -Sj  -1 

-1 

1 

-1 

__ 

(b)  Nj  -S^ 

0 

j 

-1 

0 

(c)  Sj  .  -Nj  -Nj  Sj 

0 

~+i 

14-2  Pivot  Row 

(d)  S3  =  N2  -s^  -1  S3 

-1 

1 

0 

» 

-1 

^2  Pivot  Column 
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Equation  (64a)  can  again  be  selected  as  our  cost  function. 

This  time  will  be  brought  into  the  basis  and  Sg  removed.  Again  there  cam 
be  no  increase  in  the  cost  function  because  the  minimum  ratio  is  zero. 

Once  more  expressing  the  basic  variables  in  terms  of  the  non-basic  variables 
and  transforming  the  tableau,  we  have: 


!  h 

1 

'  Si 

^2 

(a) 

s  = 

0 

-^2 

-^2 

s 

0 

.  -1 
i 

0 

■  1 

1 

(b) 

N  = 

0 

-^2 

(65) 

N 

0 

I  0 

0 

■  1 

(c) 

II 

1^5 

-^2 

-^2 

0 

-1 

!  1 

1 

(d) 

S,  = 

N, 

-1 

s. 

;  -1 

1 

Po"' 

-1 

The  first  row  of  the  tableau  has  again  lost  its  identity  as  a  cost  row.  This  is 
so  because  N2  has  a  negative  coefficient  or  cost  factor.  As  before,  use  is  made 
of  the  relationship  N2  =  -N2  and  the  necessary  changes  give: 


(a) 

S  = 
0 

N2  -S 

(b) 

N  = 
0 

N2  -S 

(c) 

Ni  = 

N2  -s 

(d) 

S3  “ 

-^2  -Si 

(66) 


s 

0 

N 

0 

h 

-1 

0 

0 

0 

1 

1 

N2  |pivot  col. 

^cost  function 
%)ivot  row 

-1 

Nl 

0 

-1 

1 

-1 

S3 

-1 

1 

0 

1 

This  time  row  one  will  be  both  cost  row  and  pivot  row  since  the  minimum  posi¬ 
tive  ratio  occurs  in  this  row.  In  this  transformation,  S  leaves  the  basis  and 

o 

therefore  assumes  the  value  of  zero. 
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Previously,  it  was  stated  that  in  maximizing  the  EN^  would  assume  the 

minimum  value  possible.  It  will  be  shown  that  the  cost  slack  can  never  assume 

* 

a  value  greater  than  zero.  Therefore,  EN^  cannot  increase  by  more  than  one 
since  must  originally  be  equal  to  minus  one. 

First  take  the  situation  where  the  cost  row  becomes  the  pivot  row.  In 
this  case,  the  cost  slack  will  leave  the  basis  and  have  a  value  of  zero. 

In  the  only  other  situation,  the  cost  row  is  not  also  the  pivot  row.  In  this 
case  if  the  cost  is  to  change,  then  the  constant  column  element  of  the  pivot  row 
must  be  non-zero.  From  (67)  it  is  seen  that  in  order  for  to  be  picked  as  the 
pivot  row 


K.  K 

#  <  ^  where  K.. 
1  c 

From  the  transformation  rules. 
Where 


E.,  K  ,  E  are  positive  numbers, 
i  c  c 

the  New  Cost  (K  )  is  received, 
c 


K  =  -old  K  ■  •'l' 

'  '  - - 


E  •  K. 

-old  K  +  -A= - - 

c  E. 


E  •  K. 

but  old  K  <  — ^ 
c  E. 


keeping  the  cost  negative. 


Concluding,  it  is  evident  that  the  cost  must  remain  negative  unless  the 
cost  row  is  picked  as  the  pivot  row,  in  ^lich  case  the  cost  goes  to  zero. 
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cost  row 


(67) 


V.  K.  i  E.  pivot  row 

11  1 


Vari.  _ 

Col. 

Con. 

Col. 


Returning  to  the  computation  and  taking  the  next  step- -once  again  the 
equations  are  put  into  canonical  form  with  respect  to  Ng  which  was  chosen  to 
enter  the  basis,  and  the  tableau  is  transformed  giving: 


1 

1 

_ i 

h 

!«! 

i  ..  i 

^2 

I  o 

(a) 

^2 

=  S 

o 

^^2 

+  1 

^2i 

* 

1  1 

i  0  1 

1-1 

i  -1 

(b) 

N 

o 

=  S 

o 

+  1 

N  : 
o  ' 

1  ; 

0 

1 

1 

1  0 

1 

:  -1 

(c) 

=  S 

o 

+  1 

X 

1  , 

-1  i 

1 

1  0 

(d) 

=»  -S 

o 

-^1  -^2 

-2 

So  i 

1 

-2 

1 

1 

1 

i  1 

Looking  at  the  equations  and  tableau,  it  is  obvious  that  a  new  cost  function 
can't  be  picked.  Observing  also  that  with  a  value  -2  remains,  which  indicates 
a  multi-core  solution  with  the  first  core  satisfying  only  the  fiist  three  states  of 
the  truth  function,  leaving  a  residue  function  »  0001  . 
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The  solution  for  the  first  core: 


1 

1  (69) 

1 

DON'T-CARE  CONDITION 

The  "don't  care"  condition  is  readily  handled  by  the  Simplex  algorithm.  This 
is  accomplished  by  removing  the  constraints  which  involve  input  conditions  which 
never  occur  or  are  of  no  consequence.  In  many  cases  tested,  this  results  in  a 
marked  reduction  in  hardware.  Further,  since  a  given  function  is  mechanized 
by  OR-ing  together  one-core  functions,  a  "don't  care"  list  may  be  generated  for 
any  problem.  This  is  accomplished  by  using  the  positions  of  the  "ones"  in  the 
sub-function,  K.  as  the  "don't  care"  conditions  in  sub-function,  K  +  1,  and  the 
positions  of  the  ones  in  sub-functions  K  and  K+  1  and  "don't  cares"  for  sub¬ 
function  K  2,  etc.  This  procedure  also  in  many  cases  resulted  in  a  reduction 
of  hardware.  For  example,  a  nine-variable  function  processed  in  the  design  of 
the  decimal  adder  was  reduced  from  a  24-core  to  a  16-core  design. 


N  = 
o 
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SECTION  m 


CONCLUSIONS 

The  study  program,  magnetic  logical  transducers,  is  concerned  with  the 
problem  of  achieving  designs  for  single  stage,  combinatorial  embodiments  of 
arbitrary  switching  functions  of  n  variables  by  means  of  magnetic  phenomena. 
Because  of  their  inherent  simplicity  and  reliability,  and  because  the  types  of 
functions  which  they  can  embody  do  not  depend  upon  physical  structure  as  in 
the  case  of  multipath  magnetic  devices,  the  study  is  confined  to  square-loop 
magnetic  toroids;  the  required  versatility  in  function  embodiment  is  attained 
by  applying  the,  at  first,  rather  nebulous  concept  of  threshold  device  logic. 

The  study  necessarily  divided  itself  into  two  complementary  phases,  the 
first  concerned  with  an  immediate  solution  to  the  design  problem  and  the  second 
with  a  theoretical  study  of  threshold  devices  having  linear  summing  input 
circuits.  The  result  of  the  first  phase  is  an  algorithm  in  the  form  of  a  com¬ 
puter  program  written  for  the  Burroughs  220  computer.  Functions  of  as 
many  as  nine  variables  are  routinely  designed  in  a  form  which  is  approximately 
minimal,  at  least  in  core  count,  by  the  use  of  this  algorithm.  The  result  of 
the  second  phase  is  a  greater  understanding  o^  the  operation  of  the  input  circuit 
of  a  threshold  device,  logic  element,  and  a  new  tool,  the  function  profile,  which 
makes  possible  the  design  of  functions  of  any  number  of  variables,  limited  only 
by  the  amount  of  labor  involved.  Furthermore,  an  approach  to  the  formulation 
of  a  logic  of  functions  realizable  in  a  single  threshold  element  is  made  possible 
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by  the  discovery  of  minimum  sets  of  circuits  which  map  all  realizable  functions 
of  a  given  number  of  variables  and  of  automatic,  means  for  describing  all  re¬ 
quired  function  transformations.  The  concept  of  fimction  modification  by 
profile  shifting  appears  to  have  application  to  the  field  of  Bionics. 


-78- 


SECTION  IV 


RECOMMENDATIONS 

It  is  recommended  that  study  continue  in  the  areas  described  below. 

THEORETICAL  STUDIES 

The  Profile  Technique  is  a  powerful  tool  which  is  used  to  study  the  input 
composite  of  a  magnetic  core  or  other  threshold  devices.  The  areas  of  interest 
include  (a)  the  formulation  of  a  realizability  test  algorithm  which  is  necessary 
and  sufficient  for  any  number  of  variables,  (b)  proof  of  all  theorems  used  in 
the  partitioning  algorithm,  (c)  the  development  of  an  algorithm  for  reducing  an 
arbitrary  solution  to  minimality,  (d)  study  of  sets  of  building  blocks  based  on 
profile  shifting. 


SIMPLEX 

Minimality 

Though  satisfactory  designs  are  achieved  by  Simplex,  they  are  not,  in 
general,  minimal,  where  minimal  is  defined  as  the  smallest  number  of  threshold 
elements  required  to  mechanize  a  given  function.  The  difficulty  is  that  the 
manner  in  which  a  switching  function  may  be  partitioned  into  threshold  functions 
is  not  unique.  It  has  been  demonstrated  that  the  threshold  functions  selected  to 
form  the  specified  switching  function  in  many  cases  depends  on  the  starting  point 
chosen  in  a  given  problem.  It  has  also  been  demonstrated  that  a  reduced  thres¬ 
hold  element  count  can  result  by  introducing  an  increasing  "don't  care"  list  to 
each  problem,  i.  e.  the  bit  positions  of  partition  k  become  the  "don't  care" 
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conditions  for  partition  k+  1.  This  is  permissable  because  the  partitioned 
functions  are  inclusively  ORed.  This  procedure  is  equivalent  to  mapping  the 
required  function  by  sets  which  are  permitted  to  overlap.  However’,  to  produce 
minimal  solutions  requires  the  development  of  criteria  which  lead  to  a  minimal 
number  of  sets  to  map  the  function.  Preliminary  thinking  indicates  that  the 
realizability  tests  based  on  the  Hamming  distance  between  bits  may  provide  a 
good  starting  point  for  such  an  analysis.  Study  is  proposed  in  this  area. 

Input  Limitations 

From  an  engineering  point  of  view,  based  on  tolerance  requirements,  an 
upper  limit  may  be  placed  upon  the  total  input  weights  and  the  threshold  value. 

The  nature  of  the  Simplex  algorithm  is  such  that  these  additional  constraints  may 
be  introduced  into  the  existing  program.  This  is  possible  since  in  forming  a 
partition  the  input  solution  is  available  for  each  intermediate  step.  A  solution 
may  be  chosen  which  contains  the  maximum  number  of  bits  and  which  is  still 
within  the  given  threshold  value  and  total  input  weight  specification.  Though  the 
overall  element  count  may  be  increased,  feasible  circuits  are  designed.  It  is 
proposed  that  the  above  feature  be  considered  for  future  study. 

Experimental  Studies 

Any  foreseeable  system  using  magnetic  threshold  logic  is  likely  to  be  com¬ 
posed  of  multi-level  nets.  For  this  reason,  it  is  necessary  to  study  the  terminal 
properties  of  a  magnetic  core  with  various  kinds  of  loads  to  develop  systems  with 
fan-in  and  fan-out  capabilities.  Newer  magnetic  materials  with  smaller  H^  and/or 
sharper  thresholds  may  provide  break-through. 

♦ 

Chao,  J.  C.,  On  Driver  Tolerances  for  Linear  Input  Logic  with  Magnetic  Cores, 
Burroughs  Corp.  internal  publication,  RM  61-25  (in  preparation). 
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GLOSSARY 


The  applied  magnetic  field,  or  magnetizing  force.  =  NlZ-t.  in 
M.  K.S.  unit. 

The  unit  of  applied  magnetic  field.  =  M  where  M  take  positive 
and  negative  integer  values  and  zero.  the  coercive  force, 

and  depends  on  the  desired  speed  of  switching. 

The  amplitude  of  the  applied  constant  current  pulse  taken  so  that  for 

N=  1,  m/I  =  1/1  =  H  . 

o 

A  matrix  operator  which  applied  to  the  solution  vector  transforms 
it  in  accordance  with  a  given  set  of  symmetry  operations  on  the  inputs. 

A  permutation  matrix  which  transforms  a  function  in  accordance  with 
a  set  of  symmetry  operations  on  the  inputs  to  the  switching  circuit 
generating  the  function. 

The  transpose  of  the  truth  table  with  an  added  column  of  "ones" 
considered  as  a  matrix. 

A  column  vector  representing  the  solution,  or  design,  of  a  switching 
circuit.  The  terms  q,  N^,  0  <  i  £  n  are  the  input  values  (weights)  of 
the  drivers,  x^. 

The  partitions  of  a  function.  A,  which  is  not  realizable  in  a  single 


The  input  variables  for  1  £  i  £  n.  For  n  »  0,  x  represents  the  bias 

o 

which  is  always  present. 

For  0  s  r  s  2^  -  1,  the  truth  value  of  the  rth  state  of  the  truth 
function. 

The  number  of  units  of  H  for  the  rth  state  of  the  function.  M  is  the 

o  r 

sum  of  the  input  values  (weights)  of  the  drivers  which  are  "on"  in  the 
rth  state. 

A  number  which  arises  in  testing  whether  a  hypothetical  profile  is 
consistent.  If  6^  2  2,  the  profile  is  not  consistent. 


